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Systematic approach to thermal leptogenesis 
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In this work we study thermal leptogenesis using non-equilibrium quantum field theory. Starting 
from fundamental equations for correlators of the quantum fields we describe the steps necessary 
to obtain quantum kinetic equations for quasiparticles. These can easily be compared to conven- 
tional results and overcome conceptional problems inherent in the canonical approach. Beyond 
CP-violating decays we include also those scattering processes which are tightly related to the de- 
cays in a consistent approximation of fourth order in the Yukawa couplings. It is demonstrated 
explicitly how the S-matrix elements for the scattering processes in the conventional approach are 
related to two- and three-loop contributions to the effective action. We derive effective decay and 
scattering amplitudes taking medium corrections and thermal masses into account. In this context 
we also investigate CP-violating Higgs decay within the same formalism. From the kinetic equations 
we derive rate equations for the lepton asymmetry improved in that they include quantum-statistical 
effects and medium corrections to the quasiparticle properties. 
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I. INTRODUCTION 

If one combines today's Standard Model of particle 
physics (SM) and that of cosmology, one finds inevitably 
that particles and their antiparticles annihilate at a very 
early moment in the evolution of the universe, leaving 
just radiation behind. The absence of a sizable matter- 
antimatter asymmetry at this epoch would imply that the 
universe as we know it could never form. The question 
about the origin of the observed asymmetry therefore 
represents a major challenge for modern physics. 

In the SM baryon and lepton number are (acciden- 
tal) global symmetries. If baryon number was also con- 
served in the early Universe a dynamical emergence of the 
asymmetry would have been impossible. In grand-unified 
extensions (GUTs) of the SM baryon number (and also 
lepton number) is explicitly broken. According to past 
reasoning, this could provide a solution to the appar- 
ent discrepancy. In the class of 'GUT-baryogencsis' sce- 
narios the matter-antimatter imbalance is generated by 
asymmetric decays of new super-heavy bosons. Anoma- 
lous electroweak processes 0, Q (sphalerons) which vio- 
late baryon and lepton number but conserve their dif- 
ference essentially eliminated the prospects for GUT- 
baryogenesis @. At the same time, it inspired the now 
widely appreciated scenarios of 'electroweak baryogen- 
esis' pi, H| and 'baryogenesis via leptogenesis' Q. Ac- 
cording to the latter scenario, the asymmetry is initially 



generated in the leptonic sector by the decay of heavy 
Majorana neutrinos at an energy scale far above the elec- 
troweak scale. Subsequently it is converted into the ob- 
served baryon asymmetry by sphalerons. The mass scale 
of the heavy Majorana neutrinos required for leptogene- 
sis [f| 0] fits together very well with the mass-differences 
inferred from observations of solar-, atmospheric- and 
reactor-neutrino oscillations. 

We focus here on the conventional, but most popular, 
high-energy (type-I) seesaw extension: 
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where Ni = Nf are the heavy Majorana fields, l a are 
the lepton doublets, (f> = ia2<j)* is the conjugate of the 
Higgs doublet, and h are the corresponding Yukawa cou- 
plings. The Majorana mass term violates lepton number 
and the Yukawa couplings can violate CP. Therefore the 
model fulfills essential requirements for baryogenesis 
They can also be realized for more complicated SM ex- 
tensions and a wide range of values for couplings and 
neutrino masses In general the right-handed neu- 

trinos do not necessarily get into thermal equilibrium and 
CP- violating oscillations between them can contribute to 
the asymmetry. This effect of leptogenesis through neu- 
trino oscillations (l3l | is crucial for neutrino-minimal ex- 
tensions of the SM (i/MSM) [H[ and poses interesting 
questions for non-equilibrium quantum field theory [l5l - 
Uj\ . In the considered scenario of thermal leptogenesis 
the heavy Majorana neutrinos experience only a moder- 
ate deviation from thermal equilibrium at the time when 
the bulk of the asymmetry is produced. Also, for a hier- 
archical mass spectrum, effects related to oscillations are 
negligible. 
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The amount of the generated asymmetry is determined 
by the out of equilibrium evolution of the heavy Majo- 
rana neutrinos. Therefore, statistical equations for the 
abundance of the neutrinos and the generated asymme- 
try are needed. The conventional approach here follows 
the lines developed for GUT-baryogenesis [l^- The CP- 
violating amplitudes for the decay and scattering pro- 
cesses involving the heavy Majorana neutrinos are com- 
puted in terms of Feynman graphs at lowest loop order. 
They are used to build generalized Boltzmann collision 
terms for these processes. Each of them contributes to 
the evolution of the distributions of Majorana neutrinos 
and leptons or, upon momentum integration, their entire 
abundances. 

However this approach is plagued by the so-called 
double-counting problem which manifests itself in the 
generation of a non-vanishing asymmetry even in ther- 
mal equilibrium. This technical issue is expression of 
the fact that the 'naive' generalization of the collision 
terms is quantitatively inexact, and inconsistent in the 
presence of CP-violation. After a real intermediate state 
(or RIS) subtraction procedure and a number of approx- 
imations, it can be made consistent with fundamental 
requirements. Nevertheless this pragmatic solution re- 
mains unsatisfactory. The requirement of unitarity guar- 
antees a consistent approximation for the amplitudes, re- 
alized by the RIS subtraction, if the statistical system 
is in thermal equilibrium. However, the deviation from 
equilibrium is a fundamental requirement for lcptogcnc- 
sis and it is not obvious how the equations have to be 
generalized for a system out of equilibrium. 

Furthermore, the CP-violation arises from one-loop 
contributions due to the exchange of virtual quanta. As 
such they seem to be beyond a Boltzmann approxima- 
tion. But the relevant imaginary part is due to inter- 
mediate states in which at least some of the particles 
are on-shell. These can also be absorbed or emitted by 
the medium and it is not obvious how such contributions 
enter the amplitudes. It is, however, clear that the in- 
fluence of medium effects on the one-loop contributions 
enters directly the CP-violating parameter and therefore 
the source for the lepton asymmetry. Their size can be 
of the same order as that of the vacuum contributions. 

Those questions can be addressed within a first- 
principle approach based on non-equilibrium quantum 
field theory (NEQFT). Several aspects of leptogenesis 
have already been investigated within this approach [T|| . 
The influence of medium effects on the generation of the 
asymmetry has been studied e.g. in |l7t [20l - [25l |. and an 
analysis with special emphasis on off-shell effects was per- 
formed in (2(| [2?} ■ The role of flavor effects as well as the 
range of applicability of the conventional approach to the 
analysis of flavored leptogenesis has been investigated in 
[28| . The resonant enhancement of the lepton asymme- 
try has been addressed within a first-principle approach 
in [r| [III Hi| [3^ . In addition, steps towards a consistent 
inclusion of gauge interactions have been taken 



In this work we use the 2PI-formalism of NEQFT to 
derive Boltzmann-likc quantum kinetic equations for the 
lepton asymmetry. In particular, we show how two-body 
scattering processes that violate lepton number by two 
units and contribute to the washout of the asymmetry 
emerge within the 2PI-formalism. This approach treats 
quantum field theory and the out of equilibrium evolution 
on an equal footing and allows to overcome the concep- 
tional difficulties inherent in the conventional approach. 
It allows us to obtain quantum-generalized Boltzmann 
equations which include medium effects and which are 
free of the double-counting problem. In other words, the 
structure of the obtained quantum kinetic equations au- 
tomatically ensures that the asymmetry vanishes in ther- 
mal equilibrium and no need for RIS subtraction arises. 
The resulting equation for the lepton asymmetry Yl is 
given by 
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Together with the 'effective amplitudes' 3 this is the main 
result of this paper. In Eq. (TT]) we introduced 



-r-PaPt.^piPj. 
ab..<^ij.. 



with £ a = 

^PaPb^-piPj 



[2ir) 4 6(p a +p b + .. - pi -pj - ..) 

x jj r A\ c/„u\ eh)- 

-fah»0--?fi)0--?f j)~]> ( 2 ) 

-(— )1 for fcrmions (bosons). Note that 
J^ab'^ij^ 3 " vanishes in equilibrium due to detailed bal- 
ance. This ensures that the asymmetry vanishes in ther- 
mal equilibrium as mentioned before. The effective am- 
plitudes contain medium effects ignored in the corre- 
sponding canonical expressions. 

We find that, in the amplitudes of the scattering pro- 
cesses medium effects are sub-dominant and can be ne- 
glected. The total decay amplitude of the Majorana neu- 
trino is barely affected as well. However, at high tem- 
peratures the available phase space shrinks when tak- 
ing gauge interactions in the form of effective thermal 
masses of Higgs and leptons into account. This leads to 
a suppression of the decay and scattering rates. Since 
the CP-violation appears as loop effect it is more sen- 
sitive to influences of the surrounding medium. Even 
though there is a partial cancellation of the fermionic and 
bosonic contributions, the CP-violating parameter is en- 
hanced by medium effects. However, the thermal masses 
reduce the enhancement and turn it into suppression at 
high temperatures. 
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We review the conventional approach to leptogencsis 
based on RIS subtraction in Sec. [II] In section Mil we 
demonstrate explicitly that in thermal equilibrium the 
success of this procedure is guaranteed by the require- 
ment of unitarity. In Sec. II VI we review the derivation of 
rate equations for total abundances and discuss in how far 
quantum statistical and medium corrections can be incor- 
porated in the reaction densities. In Scc.[V]wc review the 
application of the 2PI approach of NEQFT to leptogenc- 
sis. Equation (TT|) and explicit expressions for the effective 
in-medium decay and scattering amplitudes are derived 
within this framework in Sec. I VII We compare the results 
obtained within the 2PI-formalism to those of a conven- 
tional analysis with manual RIS subtraction. In Scc. lVIII 
we derive rate equations and the CP-violating amplitudes 
for Higgs decay within the framework of NEQFT. Finally, 
we summarize the results and present our conclusions in 

se dvnn 



II. CONVENTIONAL APPROACH 

The amount of produced asymmetry depends on the 
details of the non-equilibrium evolution of the Ma- 
jorana neutrinos as well as on the strength of CP- 
violation. The latter is usually quantified by CP- 
violating parameters [9W12I] : 



r 



where Tn^i^ and Tff are the vacuum decay rates 
to a particle or anti-particle pair respectively.For a hi- 
erarchical mass spectrum can be computed perturba- 
tively as the interference of the tree-level, one-loop vertex 




(c) 



FIG. 1: Tree-level, one-loop self-energy and one-loop vertex 
contributions to the decay of the heavy Majorana neutrino. 

and one-loop self-energy [38l - l40j amplitudes in Fig.[T] 
The contribution of the loop diagrams can be accounted 
for by effective Yukawa couplings fill ]: 



^+,ai — h>ai iH a j{h fij 7 

A_, a i = h ai — ih a j(h) h)ji f. 
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where the loop function /y is defined as 
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The first term in Eq. (|U) is related to the self-energy and 
the second term to the vertex contribution. The decay 
widths are proportional to the absolute values of the ef- 
fective couplings, rjv;->£0 = <?uj(At X + )uMi/ (32tt) and 



g w (\_\_)iiMi/(32ir) respectively, where we 
have summed over flavors of the leptons and SU {2)l in- 
dices (hence the factor g w = 2) in the final state. Since 
the phase space for the decay into particles and antiparti- 
cles is the same, one gets for the CP- violating parameter: 



x Vii , J ' + * • 



(5) 



Let us note in passing that the divergence of the loop 
function for j = i is not physical and can be removed 
by a resummation of the self-energy contribution 0- 
|42j . Here we work in a regime where the mass splittings 
| Mi — Mj\ are large enough to render effects related to 
the enhancement of the self-energy contribution irrele- 
vant (non-resonant leptogencsis). We do not require a 
strictly hierarchical mass-spectrum, however. 

To describe the statistical evolution of the lepton asym- 
metry one usually employs generalized Boltzmann equa- 
tions for the one-particle distribution functions of the 
different species fl8l . l43l . . Taking into account decay 
and inverse decay processes one writes for the distribu- 
tion function of the leptons (for a single flavor): 

p"*V/ = 5E / dni<mf{2vY5{p + k- q ) 



x [l-^l 2 Wl _^(W f )(i + /0)iV, 



(6) 



where dli^ = d 3 p/[(2ir) 3 2E p ] is the invariant phase space 
element, Sj denotes spin degrees of freedom of iVj, and T>^ 
is the covariant derivative. The corresponding equation 
for antilcptons may be obtained by interchanging I <H- I 
and (/> o (/>. CP T-invariance implies that \M\ y! .^ t ^ = 

and l-^lwi^ = \ M \ 2 i4,^N,- Furthermore, in 
thermal equilibrium detailed balance requires that (1 — 
//*)(! ~ fl q )fN\ = f?f?0- -f% t ). Subtracting the 
two relations we find for the contribution of the (inverse) 
decay terms: 

p^fe -h)= 2 • \t I du t d K i (^(p + k ?) 

i,8i 

x ail - f? )(1 + f; q )[\M\%^ erp - \M\ 2 Ni ^ $ ] . (7) 

If the decay amplitudes in square brackets differ, the 
right-hand side of Eq. ((7|) represents the (non-zero) CP- 
violating source term for the asymmetry generation. The 
total asymmetry is given by the sum over all flavors and 
SU(2)l components: = J2 a a {ne ~ ng). Neglecting 
the quantum-statistical terms, (1 — + fj,) ~ li an d 
integrating Eq. ([7]) over the lepton phase space wc obtain 
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for its time derivative: 



^0, (8) 



where K\ is the modified Bessel function of the second 
kind, Ti = r^-^ + T is the total tree-level decay 
width of Ni, and the factor gjy = 2 emerges from the sum 
over the Major ana spin degrees of freedom in ([7]), see Ap- 
pendix^ for more details. This implies that the source- 
term for the lepton asymmetry differs from zero even in 
equilibrium. On the other hand, combined with time 
translational invariance of an equilibrium state, CPT- 
invariance requires the asymmetry to vanish in thermal 
equilibrium. Thus, we arrive at an apparent contradic- 
tion. 

The generation of an asymmetry in equilibrium within 
the 5-matrix formalism is a manifestation of the so-called 
double-counting problem. In vacuum an inverse decay 
immediately followed by a decay is equivalent to a scat- 
tering process where the intermediate particle is on the 
mass shell (real intermediate state or RIS). Thus, the 



the total asymmetry til the summation over spin, flavor 
and SU(2)l components reduces to summation of the 
corresponding decay and scattering amplitudes. We will 
denote these sums over internal degrees of freedom by 3 
and call them effective amplitudes in the following. For 
the effective amplitude of £<f> — > £<p scattering one obtains 

52 



x [2{XlX + )% p:{s)P {s) + 2(/itfc)» P?(t)PiW 
+ {X\h)% P;{s)P 3 {t) + Pt(f)Pi{a)] , 

where pi 2 are the momenta of initial and final leptons 
respectively, and s and t are the usual Mandelstam vari- 
ables. The amplitude is obtained by interchang- 
ing A + A_. Note that the loop corrections to the 
Yukawas vanish for negative momentum transfer, i.e. in 
the f-channel. For this reason the above scattering ampli- 
tude contains combinations of the Yukawa couplings and 
their one-loop corrected counterparts. The propagators 
Pi are given by 



y 



(a) 



(b) 



FIG. 2: Two-body scattering process l(f> o £<f>. Both graphs 
contribute with all Ni as intermediate states. Note that we 
read (b) as t-channel contribution. 

same contribution is taken into account twice: once by 
the amplitude for (inverse) decay processes, and once by 
that for the £<j> -H- l(j> scattering processes, see Fig.|5J(a). 

Let us convince ourselves that this is indeed the case. 
Including scattering processes we have for the distribu- 
tion function of the leptons: 



x [\M\i^f^(i-mi+u) 
-m +/*)]> 



(9) 



TV) 



M 2 



i% 2 )M<r< 



(10) 



where 9 is the Hcavisidc step function. The RIS contri- 
bution appears for the flavor diagonal (i = j) terms in 
the product of the s-channel amplitudes since only in this 
case the s — M 2 terms vanish simultaneously in both Pi 
and Pi. In other words: 



8(A^A H 



(s - Mfy + (M^y 



and a similar result for 5- 



Using the defi- 



nitions of the effective couplings ([3]) and the expres- 
sion for the CP-violating parameter §5§ we find that 
(X' + X + ) 2 « (h^h)^ (1 + ti) 2 - Furthermore, for a small 
decay width, we can approximate the Brcit-Wigner prop- 
agator by a delta-function using, 



lim 



2e 



+0+ u z + 1 1 



lim 



4e 3 



0+ [lo 2 + e 2 ] 5 



2wS(uj), (11) 



where uj = s — M 2 and ti = MiTi in the considered case. 
The RIS contribution to the scattering amplitude then 
takes the form 



where the dots denote the contribution of the (inverse) 
decay processes, the sum is over flavors and SU (2)l com- 
ponents of the antilcptons and we have introduced 



, T fPaPb---PiPj- 

all ab...ij... 



Pb 



. dm 



to shorten the notation. In the unflavored regime, to 
which we restrict our analysis, the distribution functions 
of leptons of all flavors are equal. If the Majorana neu- 
trinos are close to equilibrium the difference between the 
distribution functions of the two spin degrees of freedom 
can be neglected as well. Therefore, in the expression for 



where S 



Ni->l</> 



2(piPa) 
M? 



it5 {s-M 2 ), 
MiTi 



(12) 



9w(X+X + )u2pq w g w (X\_X + ) ti M 2 



the decay amplitude squared summed over all internal 
degrees of freedom (and a similar expression for anti- 
particles). Just as one would expect, it is proportional to 
the product of the corresponding inverse decay and de- 
cay amplitudes. The additional momentum dependence 



5 



(momenta of the leptons) arises because the initial and 
final states contain fermions. Close to thermal equilib- 
rium f- e « f t fn fl q and w / w Neglecting the 
quantum-statistical terms we can write the RIS contri- 
bution to the source-term as: 

cW » 2 1 dH^^) 4 ^ fc 2 ) 

x fe 9 f</> 9 [^ij-tNi-ttj) ~ ^^iVi-^l^] ■ C-*-^) 

Taking into account that with Maxwell-Boltzmann dis- 
tributions ' f^ q = /Jy 9 in the presence of the Dirac- 
delta and performing the phase space integration using 
Eq. (|A6[) we obtain a result identical to Eq. © . 

To correct the double-counting in equilibrium we may 
therefore subtract the RIS contribution from the scatter- 
ing amplitude: 



It is defined as 



2 \pe.<j>^<t> + a e^t<^\ 

[2V ij (s)+2P*(t)P j (t) 
P*{s)P j (t)+P*(t)P j (s)) 



(16) 



Since it contains only the real part of (h'h)^j this process 
is CP-conserving. A further important washout process 





(a) (b) 
FIG. 3: Two-body scattering process U 



and similarly for the conjugate process £(f> — > £(f>. At first 
sight it might seem that the RIS subtracted scattering 
amplitudes 3 do not contribute to the generation of the 
lepton asymmetry in equilibrium, 

d t n L » 2 J efflf^ 2 * 2 (2^) 4 <5( Pl +kt-p 2 - k 2 ) 

~ "1$^>14> + a etj>-yNi-ti4 ] ' 

but also cannot compensate the asymmetry generated in 
equilibrium by the decay processes, see Eq. ([5]). However, 
upon phase space integration the difference of the unsub- 
tracted scattering amplitudes vanishes at leading order in 
h. The remaining difference of the RIS-amplitudes pre- 
cisely compensates the contribution of the (inverse) decay 
processes (©• 

The RIS subtracted scattering amplitude can be con- 
veniently rewritten in terms of a 'RIS subtracted propa- 
gator' Vij. Motivated by Eq. (fTTj) we define its diagonal 
components such, that they vanish upon integration over 
s in the vicinity of the mass pole: 



is 



[(a - Mf f + (Af^) 2 ] 2 



(14) 



Since the second of the expressions (|TT|) approaches the 
delta-function faster than the first it is common to write 
Eq. (Ti"4")l in the form 



Vu{s) P*(s)P(s) 



MiT, 



6(s - Mf) . (15) 



For i 7^ j there is no need to perform the RIS subtraction 
and therefore Vij(s) = P*(s)Pj(s). In the following we 
will also need the sum of the RIS subtracted tree-level 
scattering amplitudes. It does not contribute to the gen- 
eration of the asymmetry but plays a role for its washout. 



is 



io ^ , , scattering which receives the t- and u-channel 
contributions, see Fig. [31 By analogy with Eq. f) it i 
convenient to introduce 



2 l"AA-yU ~ "AA-tltl 

2(pu*)£ lj M.M. ^hlnl 
x [2P*(t)P*(t) + 2P*(u)P*(u 
+ P*(u)P j (t)+P*(t)P j (u)} . 



(17) 



Since the intermediate Majorana neutrino cannot go on- 
shell in the t- and u-channcl, there is no need to use the 
RIS subtracted propagator in Eq. p7|) . 

Above we have briefly reviewed the canonical approach 
to the computation of the lepton asymmetry which is 
based on generalized Boltzmann equations. Boltzmann 
equations, according to conventional reasoning, describe 
scattering processes of particles which propagate freely 
over timescales large compared to the duration of indi- 
vidual interactions. This picture seems to be consistent 
with the use of 5-matrix elements which are intended 
to describe transitions between asymptotically free ini- 
tial and final states. However, in leptogenesis the crucial 
processes ( CP- violating decays) involve unstable parti- 
cles which spoils this picture. In vacuum the amplitudes 
for such processes can be computed in terms of their 
Feynman graphs. However the 'naive' way of generaliz- 
ing the Boltzmann equation by multiplying the obtained 
amplitudes by the one-particle distributions of the initial 
states and integrating over phase space leads to incon- 
sistent equations. The origin of this problem is that the 
obtained collision terms for particle decay and inverse de- 
cay in Eq. ^ miscount the rate of particle generation. In 
a short time-interval a finite number of unstable Majo- 
rana neutrinos - formed by inverse decay of particles and 
antiparticles - decays immediately back to either particles 
or antiparticles. These contributions to particle genera- 
tion are not included in Eq. ^ where the amplitudes are 
defined in terms of Feynman graphs. For leptogenesis, 
in the presence of CP-violation, it leads to inconsistent 
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equations and must be corrected. Since the missing con- 
tribution can be constructed as the rate of a two-body 
scattering process with on-shcll intermediate state this 
issue can be addressed by the RIS-subtracting procedure 
presented above. It modifies the amplitudes for two-body 
scattering in order to cure the problem which appears due 
to the collision terms for particle decay. 1 



III. RIS SUBTRACTION WITH QUANTUM 
STATISTICS 

It is well known that unitarity has important conse- 
quences for baryogenesis and leptogenesis [H, Hf| |4(| as 
it implies restrictions for the CP-violating amplitudes. 
The issue of RIS subtraction is as well tightly related to 
unitarity as has been mentioned in e.g. [41] . As noted in 
Scc.|lTl the use of 'naive' Boltzmann equations of the kind 
© for unstable particles leads to problems such as the 
spurious asymmetry generation in the presence of CP- 
violation in the decay of the heavy neutrinos. In this 
section we show explicitly that the success of the RIS 
subtraction in thermal equilibrium is guaranteed by the 
unitarity of the S'-matrix and how it can be generalized 
to include quantum-statistical terms. The approach to 
RIS subtraction differs slightly from the one discussed in 
the previous section. 

To illustrate it we work in thermal equilibrium, f^ q = 
f t q and fl q = fl q , where 



this end we multiply Eq. (|B4[) , which follows from the 
generalized optical theorem at order h 4 , by fg 9 f^ 9 and 
integrate over cHlf. dlL 1 ?. Assuming Maxwell-Boltzmann 



equilibrium distributions we may use j 9 f^ 
the presence of the energy conserving Dirac-delta on the 
right-hand-side : 



J2 J dIi% 9 N ^) A 8{ P + k- q ) 



dnZ% P2k2 (27r) 4 s( Pl + k 1 -p 2 -k 2 ) 



■^e<t>-yi<j>. 



req req 

J i J 4> 



(20) 



We see that imposing this as a condition for the scatter- 
ing amplitudes will correctly yield dn^j dt\ cq = if we 
neglect the quantum-statistical terms in (|19[) . Equation 
(|19p represents the zeroth-order term in an expansion 
about equilibrium. Using Eq. (|20p we can therefore ob- 
tain consistent equations at this order without the need 
to specify the detailed form of E^ 7 and H. , ?J . At 
higher order (for washout contributions) we also need to 
know the sum + . see Sec. lIVl We know 

from Sec. [II] that relation ([2"0")) can be satisfied by sub- 
tracting RIS contributions from the tree-level two-body 
scattering amplitudes and taking the zero width limit: 



r 9 = (e E ° /T 



e) 



(18) 



Subtracting from the Boltzmann equation © the cor- 
responding equation for antiparticles, summing over in- 
ternal degrees of freedom of the leptons and integrating 
with J d 3 p/[(2ir) 3 E p ] we obtain in thermal equilibrium: 2 



dt ^ 



dn$ N( (2v)*8(p + k-q) (19) 

i 



+ 2 / dn™ aKa (2^) 4 Sfa +k!-p2- k 2 ) 

X \plA^ld> ~ "£<*->£<a] (1 ~ /f + fd, q )fe 9 fd, 9 



We can exploit the unitarity of the S'-matrix and CPT- 
symmetry to obtain a requirement for a consistent ap- 
proximation of the decay and scattering amplitudes. To 



l^i 'lcj)—>Ni 2MT- " JV <->-^ ' ^ ' 



Tr8{s-Mi 2 )„ 
2M,T l 



"Ni-*l<t> ■ 



(21b) 



Note that, strictly speaking, the RIS terms in Eq. (TT2")) 
include 4(pip 2 )/M 2 factors. However, upon the phase 
space integration in Eq. (|2"0"1) the two expressions give 
identical results and are therefore equal in an average 
sense. It is obvious from comparison of Eqs. (TTTH) anc ^ 09 
that the above definition of the RIS subtracted scattering 
amplitudes is not sufficient to guarantee zero asymmetry 
in equilibrium if quantum-statistical terms are included. 
However this can be achieved if we replace the vacuum 
decay width in Eq. ([21]) by the thermal one 0,113: 



1 The two pictures might seem equivalent for leptogenesis, but the 
first one implies that the Boltzmann equation for Majorana neu- 
trino decay miscounts the rate as well. This is not corrected 
by the RIS subtraction of £(j> <-> £<j> processes. However the 
corresponding correction appears at order e|, which is usually 
neglected. 

2 In order to achieve exact thermal equilibrium, in this section we 
drop the SHn^ contribution which describes the dilution due to 
the expansion of the universe. 



r < = r <(«) = / dU e p dIli(2nfS(q-p-k) 

Zg N Mi J 

x [E Ni ^ + E Ni ^] {I -ft 9 + fl 9 ). (22) 

Using the identity 1 = J ds J d A qS + (q 2 — s)5(q—p— k) and 
the fact that the (inverse) decay amplitudes are related 
by CPT-symmetry we can rewrite the RIS-contribution 
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to the second term of Eq. (fTO)) in the form: 



ds 



d 4 q 



S(s-Mf) 
2M i T i 



(2tt) 3 

J du pi dui i {2^)H{ q - Pl -k 1 ) 
J dnl 2 diii 2 (27T) 4 6(q- P2 -k 2 ) 

x(i-0i + /7)[ SjV ^- 



(23) 



The integration over s is trivial. The 5+(q 2 — s) term 
ensures that after integration over dq° the intermediate 
Majorana neutrino is on-shell, q 2 = M 2 . Using fg 9 fl q = 
InS 1 ~ fe Q + f7) together with the definition we 
can rewrite the second term of Eq. (|2"3"|) as 2gfqM i T i f e ^, 
which cancels the factors coming from RIS subtraction. 
The resulting expression reads 



pkq 



(2tt) 4 8{p + k- q) 



and cancels the first term on the right-hand side of 
Eq. ([11]). Since H^_^ = 3^_^ at 0(/i 4 ) the new 
RIS subtracted source-term for the asymmetry vanishes 
in equilibrium. 

The thermal width T i defined in Eq. (J22J) would also 
be obtained if one computes it using thermal cutting 
rules instead of the optical theorem (which applies in vac- 
uum), see Appendix[B] 3 We have seen that the unitarity 
of the 5-matrix can be employed to generalize the con- 
cept of RIS subtraction to rate equations which include 
quantum-statistical factors. As wc shall see in Sec. lVIIl 
the Majorana neutrino decay is at high temperature re- 
placed by Higgs decay if the Higgs acquires a large ef- 
fective thermal mass. In this case thermal cutting rules 
enforce relations between the amplitudes which can be 
used to obtain consistent equations, analogous to the op- 
tical theorem, see AppcndixlBl 

Note again that in Eq. ([2"3"| wc had to assume that 
the Majorana neutrinos are in exact thermal equilibrium. 
For leptogenesis this is an inconsistent assumption since 
the deviation of their distribution from equilibrium real- 
izes the third Sakharov condition and drives the gener- 
ation of the asymmetry. Not surprisingly, the NEQFT 
approach leads to a (slightly) different result for the ki- 
netic equations. However the differences between the two 
approaches enter only at an order beyond the usual ap- 
proximation as we will discuss in the next section. 



3 Note in this context that the computation of the self-energy con- 
tribution to the CP-violating parameter in thermal QFT is in 
effect only a variation of this [22l . 



IV. RATE EQUATIONS 

In this section we review the derivation of rate equa- 
tions, discuss in how far quantum-statistical and medium 
corrections can be incorporated, and compare the struc- 
ture obtained when starting from the NEQFT result 
(fTj) with the conventional form. Solving a system of 
Boltzmann-like equations in general requires the use of 
numerical codes capable of treating large systems of stiff 
differential equations for the different momentum modes 
- a cumbersome task if one wants to study a wide range 
of model parameters. In the context of baryogenesis, a 
commonly employed simplification is to approximate the 
Boltzmann equations by the corresponding network of 
'rate equations' for number densities n a or abundances 
Y a = n a /s, where s is the comoving entropy density. The 
resulting equations correspond to the hydrodynamical 
limit of the Boltzmann kinetic equations, in the comov- 
ing frame of homogeneous FRW space-time. To obtain 
evolution equations for Y a in the conventional approach, 
i.e. from Eq. wc therefore integrate the corresponding 
Boltzmann equations over g a d 3 p/ [E p (2tt) 3 ] to obtain, on 
the left-hand sides: 



dn a 
dt 



3Hn a 



sHdYa 

z dz 



where we have introduced the dimensionless inverse tem- 
perature z = Mi/T and the Hubble rate % = H\t=m 1 - 
In the homogeneous and isotropic Universe the derivative 
of the quantity Yl = n^/s can be related to the diver- 
gence of the lcpton-currcnt = (ul,0) - a quantity 
which is particularly easy to access in the first-principles 
computation - by 



sUdY L 

z dz 



On the right-hand sides we get sums of integrated colli- 
sion terms representing the effect of the different inter- 
actions. We separate contributions attributed to decays 
and scattering: 



dY L 

dz 



E 



dY, 



Li 



dz 



D 



dYL 

dz 



The decay contributions dY^J dz\ D to dY^/dz are very 
similar to the decay contributions dYjy i / dz \ D to dY^ i j dz 
and we can treat them in the same way. Reordering the 
contributions to dYi, t (jv 4 )/ 'dz\ D we nn d 



sU 



dYL 



dz 



s% 



dYL 



D 



dz 



(24) 



.D, extra 



-ppk:q 
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where the upper (lower) signs and arrows correspond to 
the rate equations for L (Ni) abundance and we defined 



sU dY Li 



z dz 
x 



D, extra 



dIll% e (2n)U(q-k-p) (25) 

(^Ni^eip ~ ^0^JVi)(l — fNi)[f<t>fe + f$fi\ > 
which corresponds to Eq. @ , as well as 



sUdYi 



D, extra 



dn 



^(27r) 4 %-*-p) 



z dz 

x [(1 - /,)(! + /,) - (1 - mi + ty] ■ (26) 

We used CPT-symmetry of the amplitudes in the deriva- 
tion of Eqs. (|25p and (|26p. Later we will see that the 
second term in Eq. (|24|) appears also in the first-principle 
approach, compare Eq. ([TJ, while the terms in Eqs. (|25| 
and ([260 are absent. This motivates the separation into 
'regular' and 'extra' terms performed in Eq. ([2"4"| . For the 
contributions attributed to scattering we get: 



sH dY L 



dz 



sH dY L 



dz 



(27) 



S. extra 



Jjjkpqr /W 



kp\qr 



-pkp\qr 

4>4 '< U+*4><j> 



kp\qr 



with 
sH dY L 



dz 



(28) 



S, extra 



(2?r) <f (fc + p - g - r) - S^_^) 

x [(1 - /,)(! + + (1 - /,-)(! + ^)/ £ / ] , 



corresponding to Eq. (fT5|) . Again, Eq. ([2"B"| does not ap- 
pear in the first-principle approach. Since in equilib- 
rium the regular terms in each of Eqs. (|24|l and (|27p van- 
ish by detailed balance we retain Eq. p^]) in the sum 
of decay and scattering contributions. The latter van- 
ishes as well in equilibrium if we adopt e.g. Eq. (|2"Tj) 
with thermal width for the RIS subtracted amplitudes 
Out of equilibrium the last terms con- 
stitute a structural difference compared to the results ob- 
tained from first-principles. This difference carries over 
to the rate equations. We will therefore analyze these 
contributions separately. 

The computational advantage of rate equations over 
full Boltzmann equations is maximized by a number of 
common approximations. In particular, assuming that 
all species are close to equilibrium and that the Ma- 
jorana neutrino distribution function (t, \q\) is pro- 
portional to its equilibrium distribution for all values 
of the momentum \q\. The temperature for all kinetic 



equilibrium distributions is set to a common value T 
while finite deviations of the chemical potential with 
small fi/T are permitted. These approximations result 
for fii/T, n^/T, n^/T <C 1 in a closed network of rate 
equations for the abundances of the form (compare with 
IE 51 51): 



sH dY L 



z dz 
Y L 



1 



(29a) 



2y; 



? (1 + c^MEifrJE) + 4< 7 f > + 4< 7 £» , 



sH dY Nz 
z dz 



veq 



where we have introduced 



<h = ^ 



(29b) 

(30a) 
(30b) 



The factor eg = 9£(3)/7r 2 sa 1.1 (we neglect the thermal 
lepton masses here) relates the chemical potential of the 
leptons to their number density, 



T 



eg ■ 



2Y, 



and the coefficient c^i takes into account that in the SM 
the chemical potentials of leptons and Higgs are related 
by ^<f> = Ctpi • M with Cff,g = 4/7 through equilibrium 
gauge, Yukawa and sphaleron interactions [48l - l50| . 

Hence, the evolution of the abundances close to equi- 
librium is roughly governed by a few average quantities 
called reaction densities which describe decay and scat- 
tering processes. We will refer to (ei7^.), (7^), (7$!) 
as CP-violating decay reaction density, decay reaction 
density and washout reaction density respectively. For 
comparison with standard results we want to maintain 
the form of Eqs. ([29)) and repeat their derivation from 
Eqs. and (|27p to obtain expressions for the reaction 
densities which take the quantum statistical factors of 
the Boltzmann equation into account. This is important 
in the present context because the thermal corrections 
to the CP- violating parameter, to be derived later, are 
of a similar kind. To this end we use that the SM gauge 
and Yukawa interactions keep the Higgs and leptons very 
close to kinetic equilibrium: 



= ( e (Ep-p-a)/T 



/« = ( 



e) 



with a common temperature Ti = = T and chemical 
potentials /if = — fig, n$ = —fxx. We shall also use f^ q for 
the equilibrium distribution functions with zero chemical 
potential defined in Eq. (TT8"[) . 

Since a chemical potential with positive sign will ap- 
pear for either the Higgs or its antiparticle, we need to 
include at least the thermal mass of the Higgs to be con- 
sistent. In the dense plasma gauge-, Yukawa- and Higgs 



9 



sclf-intcractions induce a large thermal Higgs mass of 
about OAT. With fjL^/T ~ fit/T ~ e™ c < m^/T ~ 0.4, 
the Higgs cannot acquire a condensate component. It 
is then safe to use a Bose-Einstein equilibrium distribu- 
tion function to describe the distribution of the Higgs 
particles. Using that (1 - £%) = exp((E k - fi a )/T)f a 
and hence, for a general decay collision term N <-> ab 
in the presence of the energy conserving Dirac-dclta, 4 

fah = eX P (( / i a + / i fe )/T)(l-r/a)(l-C i, /6) fP/0--f%), 

we may write: 



aLL Nab^-abi=iN- r ab*+N 



J dHf k ab (2^S(q-k-p)E abn 



x(i-r/ a )(i-e b / b ) 

x(l 



'AT 



A? 



-11 



(31) 



We can now expand the exponential in square brackets 
in the small quantity (fj, a + ji^/T. If this quantity is 
tiny at all times the integral (|3ip will not change much 
if we neglect quadratic and higher order terms. 5 For the 
zeroth-order (first) term in square brackets we use the 
linear expansion (1 - £»/,)(! - £ b / b ) « (1 - t a f?){l - 
f/rHl - (^fl q + fb 9 )] of the prefactor. The 
linear order (second) term in square brackets will appear 
preceded by just the zeroth-order factor (1 — £ a /a 9 )(l — 

et^) = fab(l - Z N f%) With 



fab =(1- Ua 9 - Ub 9 ) • 



(32) 



To write the results in a compact form we introduce de- 
cay reaction densities with quantum-statistical factors in- 
cluded: 

(X>yj> t ) = J dnf^ (2n) 4 5(q - k - p)XE N J e N «J ih (33) 
and 

(X 7 ^^(X(l-r N ^§ i ) (34) 
= f dU^(2nrS( q -k-p)XE Ni (l~r N \)r N ]f i<p , 

where E N . is the total Majorana decay amplitude. Simi- 
larly we define the scattering reaction densities as 

(X^) = J dH^;(2ir)±6(k +p-q- r)XE aMj 

x(i-r/ a e9 )(i-e fe / 6 e9 )/r/f • (35) 



4 N, a and b can be any species for which the above conditions 
apply. Here we identify N = N it a e {1,1}, b G {<t>, <j>}- 

5 By inserting equilibrium distribution functions for leptons and 
Higgs in the derivation of the CP-violating parameter we will 
neglect terms of the order e" ac (fii + fi^/T as well. 



Since E ab ^^ refers here to a CP-symmetric (tree-level) 

amplitude squared we have (Xj^j) = (X^ b ) if X is sym- 
metric as well. 

With help of Eq. (|C2[) we may separate the contribu- 
tions to dY^J dz\ D into terms proportional to A/jy. = 
(f N . — /Jy 9 ), terms proportional to Af N . ■ ^f, or just pro- 
portional to ni/T (see AppcndixICl for details): 



sH dY Li 



z dz 
sH dY Li 



A/ 



f 



N: 



(36a) 



z 



dz 
sH dY Li 



1 1 t 



Af ^4 

= -£(1 



eg 

4> 



Af. 



/r)^7^>,(36b) 



req 

■' . 



Ni 



(1"/jv 4 ) 



(1-/^) 



(36c) 



In addition we get with Eq. (|C4[1 for the extra term in 
Eq.®: 



sU dY L% 



dz 



D. extra 



= 2((1 - f Ni )ea^) . 



Equations (JM]) describe the generation of a net asym- 
metry due to out of equilibrium decays of heavy Majo- 
rana neutrinos. Once ne/T has a non-zero value, there 
will be a slight difference in the decay rates to particles 
and antiparticles respectively which is not due to CP- 
violation in the decay amplitude, but due to the pres- 
ence of slightly different occupation numbers of leptons 
and Higgs in the final states of the decays. At linear order 
this combined effect of blocking and stimulated emission 



is accounted for by Eq. (|36b|) . Depending on the 'typical' 
sign of (c^ef^ 1 — f(' q ) it can add to or diminish an existing 
asymmetry. Finally, Eq. (|36cp describes washout due to 
inverse decays. In Sec.|V] we will see that the functional 
dependence on f^ = (I + f eq — f eq ) in the integrated 
collision terms is the same as that encountered in the 
CP-violating parameter ti itself. 

Considering the last two terms in Eq. (|27|) we find for 
the scattering contributions: 



sUdY L 



z dz 



s — 



-4f(l+«*)[<T#> + <-$>], (37) 



where we defined the 'RIS subtracted reaction density' 
(7^). If we adopt the amplitudes defined in Eq. (j2"T]) in 
the framework of RIS subtraction, it is given by 



^(l + e i)"iVi x , „ ,2\1 e<t> 



4E 



Ul, 



Note that the contribution proportional to ef is of higher 
order in h. Furthermore, we get for the extra term: 



sH dY L I 

z dz I S, extra 



-2E(( 1 -/a>-)^,) 



(38) 
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Here we used 

{16w) 2 pqE m S(s 



Mi 



a 2 M-E 



~7<>j, 



hi) 



We have written Eq. ([38]) schematically in order to show 
how it compares to other washout terms. Note that the 
extra terms indicate that there will be a slight difference 
between the equations obtained in the 2PI approach and 
those obtained with RIS subtraction at finite tempera- 
ture. Comparing Eqs. (|3"8"|) and ([3"T)l we see that the first 
term in the former equation will cancel the latter con- 
tribution in thermal equilibrium (f N . = ) if the de- 
cay contributions are summed up. The second term in 
Eq. (|38p is due to quantum statistics. Since it is pro- 
portional to ti/i/T it can be large only if is large (as 
in the case of resonant leptogencsis). Anticipating our 
knowledge about the structure obtained within NEQFT, 
we will ignore the extra term in what follows. 

At the time being, everything is still exact with re- 
spect to deviations of f N . from equilibrium. This dis- 
tribution is necessarily distorted due to the fact that 
it is subject to conflicting equilibrium conditions corre- 
sponding to the decay into particles and antiparticles, 
by the effects of the expansion and, possibly, due to 
non-equilibrium initial conditions. In order to obtain 
the full momentum-dependent distribution function we 
would have to solve the corresponding full kinetic equa- 
tions however [2(1 HI M^M ■ 

To proceed we shall as usual assume that the devia- 
tion of the Majorana neutrinos from equilibrium is small. 



This allows us to neglect the Af N . contribution (|36bl 



(J. The extra terms 



and to replace f N . — > fjj. in Eq. 
cancel at this level of approximation up to the quantum- 
statistical term. In order to bring the remaining source- 
term Eq. (|36a[) into the conventional form, we need to 
assume that the non-equilibrium distribution of the Ma- 
jorana neutrino is proportional to its equilibrium value 
(with momentum independent prefactor) 6 

f = Tljv ' f e i 

JNi — eq JNi ■ 
n N, 

With this approximation we can write 

'Y N , 



v&q 

Ni 



lUX-y) 



The total contribution to the evolution equations for the 
lepton asymmetry is then given by 



dY L 

dz 



E 



fdY L 



\ dz 



dY Ll 




dY L 




dz 


Of) 


dz 


s 



This amounts to the assumption that its shape can, in terms 
of its quantitative effect on the integrated collision terms, ef- 
fectively be captured by a Maxwcll-Boltzmann distribution with 
(small) 'pseudo-chemical potential' |44| . Strictly speaking, it im- 
plies that we need to revert to a classical distribution function 
for the Majorana neutrinos. 



i.e. we obtain Eq. (|29a[) . We see that, at this level of 
approximation, there are no contributions due to extra 
terms apart from those which cancel due to the RlS-sub- 
traction. Quantitative differences can arise if the devia- 
tion of the Majorana neutrinos from equilibrium is large 
or ti is of order 1. For the evolution of the Majorana 
neutrino we obtain with Appcndix[Cj similar to Eq. (|36p , 
the decay contributions 



sU dY Ni 



z dz 
sH dY N . 



z 



dz 



A/* 



A/ 



J n 



(39a) 



77r(ei(c^/ - ft ) , 
1 Jn 



q iNi) 



sU dYj 



Ni 



dz 



(1 /atJ w 

a-r N q f N * 



(39b) 
(39c) 



and for the extra term in Eq. ([24 



sH dY, 



Ni 



z dz 



D. extra 



"2 — (1 + C^)^e ^ 1 _ feq^ N t 



*%H°#f? - fn^g"i8 t ) • (40) 



Neglecting again Af N and contributions we ob- 
tain 



dY, 



N: 



dY, 



n, 



dz 



dz 



AA 



i.e. Eq. (|29b[) . If higher order contributions are taken 
into account, we get a difference between the conven- 
tional equations and those derived in the 2PI-formalism. 



Ignoring the contribution (J39bj) and the second term in 
Eq. ([4"0"[) . which are due to quantum statistics, we obtain 
a contribution 



=F 



Y L 
2Y e eq ' 



'■1 'nJ 



to dY^/dz. Here the upper sign applies if the extra 
terms are included and the lower sign if not. This can 
therefore result in the inclusion of this term with wrong 
sign even if quantum statistics are neglected, compare 
e.g. 



41 



The reaction densities for decay, (ei7^.), (7^.), (7^), 
and scattering, (7^ ), (7^), represent the hydrodynami- 
cal coefficients which govern the evolution of the number 
densities (abundances). We will compute them numeri- 
cally once the additional medium dependence of the am- 
plitudes (in particular the CP-violating parameters) has 



7 The origin of this difference is that no RIS subtraction alike is 
performed for the Boltzmann equation of the heavy Majorana 
neutrino. 
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been derived. In addition, it is useful to define a ther- 
mally averaged CP- violating parameter as 



(jag, 



which equals if it is momentum independent, such as 
in the zero temperature case, but will differ once ther- 
mal effects are included. This quantity is meaningful 
for the comparison with conventional results because it 
takes into account that the deviation of the Majorana 
neutrino abundance from equilibrium, which appears in 
the source-term for the lepton abundance, is influenced 
by the ( CP-conserving) decay reaction density in the de- 
nominator. 

Inserting conventional vacuum amplitudes in Eqs. (1291) 
with Eqs. (I33|) and (|34f) and dropping quantum-statistical 
factors one obtains the conventional results for the reac- 
tion densities. For the readers convenience we quote them 
here. For the decay reaction density we obtain 



<7#> = <7ft> 



9N 

2tt 2 



MfTiTKx 



Mi 
T 



(41) 



and (e, 7^ 



e, (7jy. ) , see Appendix[21 For the two-body 
scattering the reaction density is given by 



T 
64tt 4 



OO _ 



(s), (42) 



where a(s) is so-called reduced cross section: 



(43) 



subtracted cross section' of the £(f> o £<f> process: 



(45) 



x < x 



(x — cii)(x — cij) + (1 — 2Sij)aia,jCiCj 
[(x - ai) 2 + (flic,) 2 ] [(a; - Qj) 2 + (ajCj)" 2 



2^±* In 
dj — ai 

x — a 



x + ai 



In 



(x ~~ a^ 2 + (a,c,) 2 
x — a 4 



x — (x + aj) In 



x — (x + a^ In 



x + a,j 



X + Gj 



x + ai 



(x - a,j) 2 + {ajCj) 2 

The reduced 'cross section' f|45[) is negative® in the vicin- 
ity of the mass shells, x « aj. This is due to the — (MiTi) 2 
term in the numerator of the RIS subtracted propaga- 
tor HH). Note that because we have not approximated 
this term by the Dirac-delta the structure of Eq. (1431) is 
slightly different from the one usually used in the litera- 
ture [|l|. 

V. NON-EQUILIBRIUM QFT APPROACH 

In this section we briefly review the description of lep- 
togenesis within non-equilibrium quantum field theory 
[54M57| . This framework has been shown recently to be 
suitable for the derivation of quantum dynamic equa- 
tions for the lepton asymmetry within a first-principle ap- 
proach, and to incorporate medium, off-shell, coherence 
and pos sibly further quantum effects in a self-consistent 
way [19I4281 . [58M63} . We continue these efforts by deriv- 
ing consistent quantum corrected Boltzmann equations 
that describe the generation and washout of the lepton 
asymmetry and include the (inverse) decay as well as 
scattering processes mediated by Majorana neutrinos. 



For the it -H- <fi(j) process it reads 



In 



a, (x + aj ) 
aj(x + ai) 



2x 



■In 



(x + a t )(x + aj) 



(44) 



where we have replaced s by x = s/Al 2 and introduced 
dimcnsionless quantities aj = M 2 /M 2 and c, = Ti/Mi = 
(h'h)ii/87r. The case i = j is included in this expres- 



sion in the limiting sense aj — > ai. Note that Eq. (|44|) 
only contains the real part of (h^h) 2 . The contribution 
of the imaginary part vanishes because 3(/i^/i) 2 j is anti- 
symmetric with respect toif> j whereas the sum in the 
curly brackets is symmetric under this transformation. 
The integration of Eq. (TIB)) yields for the reduced 'RIS 



A. CTP and propagators 

The lepton asymmetry is given by the (i = 0-compone- 
nt of the expectation value of the lepton-current operator: 

i£(s) = <E£(s)V£(aO>. 

a, a 

It can be expressed in terms of the leptonic two-point 
function. We define the two-point functions for the Higgs, 



8 Note that a' is not a physical cross section but denotes the con- 
tribution to the reaction density arising from the difference of 
the full and the RIS term. We stress that all physical rates are 
manifestly positive, e.g. the washout term, to which cr' yields 
a sub-leading correction that is relatively suppressed by Yukawa 
couplings. See also [4ll . 
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lepton and Majorana fields with time arguments attached 
to the closed time path (CTP) shown in Fig. [5] by 

A o6 (x,y) = (T c 0°(x)^ b (t/)), (46a) 
S^{x,y) = {T e t a {x%{y)), (46b) 
j7 ij (x,y) = {T c N i (x)N j {y)}, (46c) 

where the sub- and superscripts refer to SU(2)l and fla- 
vor indices and Tq denotes time-ordering with respect to 
the CTP. We will frequently use matrix notation for the 

' , =^ 

tiriit tmax 

FIG. 4: Closed time path. 

flavor indices, where e.g. 5? denotes the flavor-matrix 
S^ 3 , etc. Using the definition (|46b[) we find for the 
lepton-current: 

j£(*) = -E tr [vsrM] ■ 

a .a 

Two-point functions G(x,y) (where G stands for A. S 
or ,5^) defined on the CTP can be decomposed into a 
spectral function G p and statistical propagator Gf' 

G{x,y) = G F {x,y) - |sgn c (z° - y°)G p (x,y) . (47) 

The signum function is either +1 or —1 depending on 
whether a; or y° occur 'later' on the contour C. Gf and 
G p encode information on the state and the spectrum of 
the system, respectively. For example, for the leptons 
they are given by 

s$ p (x,v) = i{[r a (x),%(v)] + )> 

where [.,.]± denote (anti-)commutators. Statistical and 
spectral functions of Majorana neutrino and Higgs can be 
expressed similarly, with + and — exchanged for bosons. 
Although there are only two independent two-point func- 
tions for each species, it is convenient to introduce addi- 
tional combinations of them, namely the Wightman func- 
tions 

G^,y)^G F (x,y)T^G p (x,y), (48) 

as well as retarded and advanced functions, 

G R (x,y) = Q(x -y°)G p (x,y), (49a) 
G A (x,y) = -@(y°-x°)G p (x,y). (49b) 

From the above definitions one can see that the differ- 
ence of the retarded and advanced propagators gives the 



spectral one, whereas the sum yields the hermitian prop- 
agator G h {x,y): 

Gr(x,v) - G A {x,y) = G p (x,y) , (50a) 
G R (x,y) + G A {x,y)=2G h (x,y). (50b) 

Finally, we will also need the CP conjugated propagators 
on the CTP: 

& ab (x,y) = A ba (y,x) , (51a) 
S a J{x,y) ^CP)S p b :{y,x) T (CP)- 1 , (51b) 
y ij (x,y) = (CP^'faxfiCP)- 1 . (51c) 

Here x — (a; , —x), C = ij 2 ^ and P = 7 are the charge 
conjugation and parity matrices, respectively, and the 
transposition refers to spinor indices. CP conjugated sta- 
tistical and spectral functions immediately follow from 
the above definition by inserting the decomposition (|47[) . 



B. Kadanoff-Baym equations for leptons 

The time-evolution of the two-point functions is de- 
scribed self-consistcntly by the Kadanoff-Baym (KB) 
equations. These equations can be obtained from a 
variational principle using the so-called 2PI effective 
action (64|. The resulting equations of motion have 
the form of Schwinger-Dyson equations for the non- 
equilibrium propagators formulated on the CTP: 

S- 1 (x,y) = So 1 (x,y)-£(x,y). (52) 

Here S~ 1 (x, y) is the inverse of the full lepton propagator 
in coordinate space, and Sq (x, y) is the inverse of the 
free lepton propagator, 

So 1 (x,y) = 6 ab S^0J c (x-y)P L . 

The information about the interaction processes is en- 
coded in the self-energies £. They can be obtained by 
cutting one line of the 2PI contributions to the effective 
actions. The two- and three-loop contributions are pre- 
sented in Fig. [51 

The KB equations can be obtained by convoluting 
the Schwinger-Dyson equation with the full propagator, 
which yields: 

$ x S^{x,y) = 

S a ^6 c (x-y) + f d 4 z^(x,z)S^(z,y). 
Jc 

Here J c d 4 z = J c dz° J d 3 z. After decomposing the re- 
sulting equation into statistical and spectral components, 
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equations ([53")) one obtains 9 : 




s 

(a) 



(b) 




(c) 



(d) 



FIG. 5: Two- and three-loop contributions to the 2PI effective 
action and the corresponding contributions to the lepton self- 
energy. Note that the propagator lines used here denote full 
resummed propagators in contrast to those employed in the 
previous Feynman graphs. The contributions (a) and (c) to 
the 2PI effective action are known as 'setting-sun' and 'mer- 
cedes' diagrams respectively. 



one obtains: 



i$ x Sf{x,y) 



' d 4 zE a /(x,z)SjP(z,y) 



d i zTfp{x,z)Sf{z,y) i (53a) 

,0 

' d 4 z^(x,z)Sf(z,y). (53b) 



The equations for Majorana and Higgs propagators have 
a similar structure, with the Klein-Gordon instead of 
the Dirac operator for the latter. The Schwinger-Dyson 
equations (f5!Zj) and the corresponding Kadanoff-Baym 
equations (|53[) are formally very similar to the Schwinger- 
Dyson equation in vacuum. However, out of equilibrium 
the propagators depend not only on the relative coor- 
dinate s — x — y, but also on the central coordinate 
X = (x + y)/2, which makes their solution much more 
involved. In contrast to the Schwinger-Dyson equation in 
vacuum, the KB equations determine the spectral prop- 
erties of the system including medium corrections, as well 
as the non-equilibrium dynamics of the statistical prop- 
agator self-consistently. Since the latter represents the 
quantum field theoretical generalization of the classical 
particle distribution functions, KB equations can be seen 
as the quantum field theoretical generalizations of Boltz- 
mann equations. 

As pointed out above, an equation of motion for the 
lepton asymmetry can be derived by considering the di- 
vergence of the lepton-current T> ll j^{x). Using the KB 



ZVJlM = ~9 W lim (2>£ + 2?£)tr [ %l S aa (x, y)] 



9w i 



ztr[Zf(x,z)S F a (z,x) 



Pot i 



T.f (x,z)S^(z,x)-Sf(x,z)^ F a (z,x) 



+ Sf(x,z)^ p a (z,x)}. 



(54) 



Here summation over repeated indices is implicitly as- 
sumed. The two equations above represent the quantum 
generalization of the Boltzmann equation for the lepton 
asymmetry. Thus, they may be considered as the mas- 
ter equations for a quantum field theoretical treatment 
of leptogenesis [Til HI] ■ 

The dependence of the two-point functions on the rel- 
ative coordinate s is characterized by the hard scales like 
the Majorana neutrino mass M\ or the temperature T of 
the surrounding plasma. In contrast to that, the varia- 
tion with the central coordinate X is given by the macro- 
scopic time-evolution of the system, e.g. the Hubble rate 
H or the Majorana decay rate T. Therefore, it is possi- 
ble to perform an expansion in slow relative to fast time- 
scales, i.e. in powers of e.g. T/M\ or H/T. Technically, 
this can be realized by a so-called gradient or derivative 
expansion with respect to X , and a Fourier transforma- 
tion with respect to s, known as Wigncr transformation, 
see Appcndix[D] for more details. Then, to leading or- 
der in the gradients, the evolution equation (JSH) for the 
lepton asymmetry becomes Markovian, and after some 
straightforward algebra, can be written as 



d 3 p 
(2^)3 



(55) 



x tr{ [E^(t,p)Si a (t t p) - ^(t,p)S^ a (t,p)] 



- [Sf (t,p)5^(t,p) - £<l a (t,p)sf(t, P )}}. 

Note that it is possible to investigate higher orders in 
the derivative expansion systematically [24| . In Eq. (|55p 
all two-point functions are evaluated in Wigner space, 
where p is the physical momentum (65^ that corresponds 
to s. For a spatially homogeneous system (like FRW) the 
two-point functions depend only on the time coordinate 
t = X , and on the momentum p, because of spatial 
translational invariance. Strictly speaking, this is true 
only in the rest frame of the medium (comoving frame). 
In a general frame the two-point functions depend on 



9 We assume here that in FRW space-time the effects of the Uni- 
verse expansion can be captured, to the required accuracy, by in- 
troducing the invariant integration measure £^ 4 z = yj—gd^z and 
using the covariant derivative X> M . As has been demonstrated in 
[65ll , this is the case for scalar fields. A manifestly covariant gen- 
eralization of center and relative coordinates X and s to curved 
space-time can be found in [66(1 . 
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X ■ u, where is the four- velocity of the medium. The 
latter satisfies the normalization condition u^u^ = 1, and 
is given by u — (1, 0, 0, 0) in the medium rest frame. 

In order to allow for a physical interpretation of 
Eq. (|55[) we have written it such that the integration 
is over positive frequencies only, and expressed the lep- 
ton propagator and self-energy in terms of the Wigner 
transformed Wightman functions Eq. (|4"5| . In ther- 
mal equilibrium, the Wightman functions depend only 
on the momentum p and satisfy the Kubo-Martin- 
Schwinger (KMS) relation Gf(p) = ±e pu / T G* < ' l (p) for 
fermions/bosons, respectively. When inserting the KMS 
relations for propagators and self-energies into Eq. (f55j) . 
one immediately finds that the divergence of the lepton- 
current vanishes in thermal equilibrium as it should (see 
also p3j). In other words, the quantum equation for 
the lepton asymmetry is in accordance with the third 
Sakharov condition. We emphasize that it is not neces- 
sary to apply RIS subtraction to obtain this result within 
the CTP approach [H HH, H| . 

The four terms on the right-hand side of Eq. ([55)) may 
be interpreted as gain and loss terms of leptons and anti- 
leptons respectively [20l l2l|. In particular, one may de- 
fine generalized lepton distribution functions /" (t, p) via 
the so-called Kadanoff-Baym ansatz 

Sf = (1 - ff)Sf , Sf = -ffSf. (56) 

Thus the contribution on the right-hand side of Eq. (|55|) 
that contains S*< Q corresponds to the lepton loss term, 
while the contribution proportional to S> represents the 
lepton gain term. Analogous definitions relate the CP 
conjugate propagators with the anti-lepton distribution. 
Note that the KMS relations ensure that in equilibrium / 
approaches the Fermi-Dirac distribution J"' 3 — > S a ^ fp D - 
The flavor off-diagonal components encode coherent fla- 
vor correlations [28| . In the unflavored regime considered 
here = S a P f e and = S afi S. In the quasiparticle 
(QP) approximation, the spectral function is given by 

Sf(t,p) = (2tt) sign(p°) 6(p 2 - ml) 8 afi P L fP R 

= S p S af3 P L pP Rl (57) 

where we assume that leptons obey conventional dis- 
persion relation and mi is the effective thermal mass. 
These assumptions might be modified in the presence of 
a medium (Tfl, H3, HI . 

Due to the presence of the Dirac-delta-function in 
Eq. (|57p the integration over p° in Eq. (j55|) is trivial and 
leaves only the integration over spatial momenta of on- 
shell leptons. Therefore the right-hand side of Eq. (f55|) 
can be interpreted as a difference of two (integrated) 
Boltzmann-likc equations - one for the particles and one 
for the antiparticles [2l| . According to the physical inter- 
pretation of Eq. (f55j) in terms of gain and loss terms, the 
Wightman components of the lepton self-energy and of 
its CP conjugate are the analogs of the collision integrals. 
Since we limit our analysis to the unflavored regime, it 



is convenient to perform the summation over the flavor 
indices: S al3 T, al3 = = S. Then the one-loop contri- 
bution, sec Fig.[5j(b), takes the form: 

Eg ) (*.f) = "/ dnidn 4 q (2n) 4 5(p + k-q) 

x {rth)ji P R &g(t, q)P L ^(t, k) , (58) 

where dll^ = d i q/{2ir) i . The explicit expression for the 
two-loop contribution is rather lengthy and it is conve- 
nient to split it into three distinct terms: 

e£> = + eL 2 ' 2) + S^ 3) . (59) 

K K < < K ' 

The first term on the right-hand side reads 

4 21) (^) = / dU t dU t ( 27r ) 4 ^ + k - ( 6 °) 
X [(rfh) in (rfh) jm A mn (t, q, k)P L Cy^{t, q)P L A^(t, k) 
+ {rfh) ni (rfh) mj P R ^(t,q)CP R V nm (t,q,p)A^(t,k)] , 

where we have introduced two functions containing loop 
corrections: 

A mn (t, q ,k) = J du 4 ki dn 4 k2 dn 4 k3 

x (27r) 4 <5(<7 + ki + k 2 ) {2n) i 5(k + k 2 - k 3 ) 
x [P R y R n (t, -k 3 )CP R &£(t, k 2 )A A (t, h) 

+ p R yp n {t, -k 3 )cp R sl(t, k 2 )A A (t, h) 
+ p R y^ n {t,-k z )CP R s T A {t,k 2 )A F {tM)} , 

and V nm (t, q, k) = P A^^t, q, k) P to shorten the nota- 
tion. Comparing Eqs. (|55|) and ([50]) we see that they have 
a very similar structure. First, the integration is over mo- 
menta of the Higgs and Majorana neutrino and the delta- 
function contains the same combination of the momenta. 
Second, both self-energies include one Wightman propa- 
gator of the Higgs field and one Wightman propagator of 
the Majorana field. Upon the use of the Kadanoff-Baym 
ansatz the Wightman propagators can be interpreted as 
cut-propagators which describe on-shell particles created 
from or absorbed by the plasma [69j . On the other hand, 
the retarded and advanced propagators can be associated 
with the off-shell intermediate states. We therefore con- 
clude that Eqs. (|55|) and (|60[) describe (inverse) decays of 
the heavy neutrino into a lepton-Higgs pair. 

The second term on the right-hand side of Eq. ([55]) 
contains two Wightman propagators of the Higgs field 
and one Wightman propagator of the lepton field. The 
Majorana propagator appears only in the intermediate 
state: 

Xf 2) (t,p) = J dU 4 p2 dUi i dUi 2 (2n) 4 S(p + k 1 -p 2 -k 2 ) 

x {rfh) ni (rfh) mj [p R y R 3 {t, P2 + k 2 )CP R s^(t, - P2 )p l 

x C<y A nn {t 1 p 2 ~k l )P L A^{tM)A^{t,~k 2 )} . (61) 
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We therefore conclude that this term describes lepton 
number violating scattering processes mediated by the 
heavy neutrino. Finally the last term in Eq. ([59|) contains 
two Wightman propagators of the Majorana field and 
one of the lepton field, whereas the Higgs field is in the 
intermediate state: 

^' 3) (t,P) = J dn P2 dU kl dH k2 (2TT) 4 6{p + qi -p2-q2) 

x (h^h) ij (h^h) lk [P R y^(t,-q 1 )CP R Sl{t,p 2 )P L C 

x ^{t,q 2 )P L A A {t,-q 2 -p 2 )A R {t,q 1 -p 2 j\ . (62) 

Therefore it can be identified with the Higgs mediated 
scattering processes. These conserve lepton number and 
do not contribute to generation of the lepton asymmetry. 

The CP conjugate of the Wigner transforms can be 
obtained using Eq. (|5ip . In practice this amounts to re- 
placing the propagators by their CP conjugate and the 
couplings by their complex conjugate in the above expres- 
sions. For instance for the CP conjugate of the one-loop 
self-energy we find: 

Zf(t,p) = - j dnt dHfan^Sik+p- q) 

x (h^h)* l P R y^(t,q)P L A § (t,k). (63) 

Expression for the CP conjugate of the two-loop lepton 
self-energy can be obtained in a similar way. 

For the Higgs propagators in the above self-energies we 
can also use the Kadanoff-Baym ansatz, 

A > = (l + /^)A p , A < =/ A p , (64) 

and the simple quasiparticle approximation for the spec- 
tral function, 

A p (t, k) = (2tt) sign(fc°) S(k 2 - ml) , (65) 

where is the effective thermal mass. Effects of the 
finite thermal Higgs mass will be studied in Sec. lVIII 

VI. MAJORANA CONTRIBUTION 

In this section we will analyze the lepton number and 
CP-violating (inverse) decay of the Majorana neutrino 
as well as the two-body scattering processes mediated by 
the heavy neutrino. In particular, we will derive expres- 
sions for the in-medium CP-violating parameters, decay 
widths and scattering amplitudes. We will also explicitly 
demonstrate that the obtained equation for the lepton 
asymmetry is free of the double-counting problem. 

A. Decay at tree-level approximation 

In the previous section we have used the Kadanoff- 
Baym ansatz and quasiparticle approximation for the 



Higgs and lepton fields. Let us now assume that similar 
approximations also hold for Majorana neutrinos. That 
is, we assume that in Eqs. (|58|) and ([63")) the spectral func- 
tion S^ %3 is diagonal in flavor space and can be approxi- 
mated by 

= (2tt) sign( g )% 2 - Mf) 5 l3 (0 + M t ) , 

and that it is related to the Wightmann components via 
the Kadanoff-Baym ansatz: 

•/'.' = (l - fn)*?, ■/" = -fN t J? ■ 

Substituting Eqs. ([58]) and ftJo) in Eq. ([So]) and making 
the above approximations we find after some algebra that 
the lepton-current can be represented in the form: 

i J 

where J- have been introduced in Eq. ^ and we have 
defined: 

%« Ni = 9w{^h)utr[{f + Mi)P L f] , (67a) 
%^ Ni = 9 w (rfh)utr[($ + Mi)P L }) . (67b) 

The superscript "I" stands for 'tree- level'. The expression 
((66)) strongly resembles the Boltzmann equation. There- 
fore the functions and Sj^^ can be interpreted 
as effective in-medium amplitudes squared, summed over 
internal degrees of freedom, for the decays into leptons 
and antilcptons respectively. The two effective ampli- 

\ \ 

FIG. 6: Tree-level contribution. 

tudes (|67]) can be replaced by the total decay amplitude 
and the CP-violating parameter. Using Eq. (|o7|) we find 
that within the used approximations the resulting decay 
amplitude coincides with the outcome of the vacuum cal- 
culation, E N , = 4g w (h^h)a(pq), and that ej = 0. 

In the presence of a nonzero lepton asymmetry fp ^ /| 

and U^f$. Therefore Ff^m ^ and this leads 

to a washout of the asymmetry. Despite the fact that 
Eq. ((66]) correctly describes the (leading-order) washout 
processes, it fails to describe processes which gener- 
ate lepton asymmetry: in the considered approximation 
e, = because the CP-violating effects, which are re- 
quired to produce the asymmetry, are of fourth order 
in the Yukawa couplings of the Majorana neutrino. In 
Eq. (|67[) we have taken into account only terms quadratic 
in the coupling. In other words, this approximation cor- 
responds to the tree-level approximation in the canonical 
approach. 
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Terms of higher order in the couplings emerge from 
three- and higher-loop contributions to the lepton self- 
energy, see Eqs. (|59|) - ([6T|) . as well as from expansion of 
the full Majorana propagators entering the self-energies. 

B. Equilibrium solution for Majorana propagator 

In order to define an effective CP-violating parameter 
and decay width that incorporate medium corrections we 
have to identify the quasiparticle excitations in the sys- 
tem. To perform this analysis we follow the discussion of 
the self-energy contribution within a toy-model as pre- 
sented in [21) . As has been demonstrated there, it is 
important to take the matrix structure of the Majorana 
propagator in flavor space into account. Our starting 
point is the Schwingcr-Dyson equation for the Majorana 
two-point function: 

y- 1 {x,y)=^ 1 {x,y)-il{x,y). 

Let us split the self-energy into diagonal and off-diagonal 
components in flavor space and introduce a diagonal 
propagator S defined by the equation: 

S- 1 {x,y)=S^{x,y)-tl d {x,y), (68) 

where S is the free propagator and IF denotes the di- 
agonal components of the self-energy. The poles of the 
diagonal propagator define the quasiparticle excitations. 
It can be shown that the dynamics of these is described 
by a Boltzmann-likc quantum kinetic equation. 

Inserting this decomposition into the Schwingcr-Dyson 
equation we find, using matrix notation: 

y-\x,y) =S- 1 (x 1 y)-fl'(x,y) 1 (69) 

where II denotes the off-diagonal components of the 
self-energy and the full neutrino propagator includ- 
ing flavor-diagonal and flavor off-diagonal contributions. 
Multiplying Eq. ([6"9"[) by &* from the left, by S from the 
right and integrating over the contour C we obtain a for- 
mal solution for the full non-equilibrium propagator: 

y(x,y) =S(x,y) 

+ J ^u@\y(x,u)i\(u,v)S(v,y). (70) 

After decomposing the propagators and self-energies into 
the spectral and statistical components, we can rewrite 
Eq. (JTUJ) in the form: 

y F{p) {x,y) = S F(p) (x,y) - J ® 4 u^ve(u )9(v ) 

x [,y R (x,u)U R (u,v)S F(p) (v,y) 

+ <?r(x, u)tl F{p) (u, v)S A (v, y) 

+ ^F{ P )( X : u )^'a( u ^ v)S a (v, y)] . (71) 



Here we are using the retarded and advanced propaga- 
tors defined by Eq. ([4"9"[) . so that the integration can be 
extended to the whole wu-plane. Using their definitions 
and Eq. (|7ip , we can also derive formal solutions for the 
retarded and advanced propagators: 

y R(A) (x, y) = S R(A) (x, y)-J @ 4 u@ 4 v9(u°)e(v ) 

X ^R(A)( X ' U )^R(A)( U ^ V )^R(A)( V ^V) ■ ( 72 ) 

Next we Wigner transform Eqs. ([TTj) and (f72"[) and per- 
form the leading order gradient expansion as has been 
outlined in Sec.[V] Combining both results, we find for 
the full statistical and spectral propagators and the cor- 
responding causal two-point functions of the system in 
to equilibrium: 

F{p) = ^F{p) ~ ^R^-R^F(p) 

~ ^R^-F(p)^A~ ^F^y^-A^Ai (73a) 

R(A) = Sr(a) - ^ r{A)Rr(A)Sr{A) > (73b) 

where all propagators and self-energies are evaluated at 
the same point (A, q) in configuration space. We can ex- 
press the full statistical and spectral propagators in terms 
of the diagonal ones and the off-diagonal self-energies, 

F(p) = ©fi [Af(p) ~ SrH-f(p)Sa\ ©a > (74) 

where Q R and Q A are defined by Q R = (l+S R fl R ) 1 and 

Q A = (1+11^5^) 1 respectively, with I being the 4nx4n 
unit matrix in the Dirac and flavor space of the n gen- 
erations. Solution Eq. ([74"[) reduces the dynamics of the 
full statistical and spectral propagators to the dynam- 
ics of two quasiparticle excitations. Their masses, decay 
widths and CP-violating parameters arc determined by 
the medium and the abundances are described by the cor- 
responding one-particle distribution functions. Strictly 
speaking, the solution (|74[) is valid only in thermal equi- 
librium. However, we assume that it also holds for small 
deviations from equilibrium. 

To consistently analyze processes of the fourth order in 
the coupling one has to use so-called extended quasipar- 
ticle approximation (eQP) for the statistical propagator 
and spectral function [70j-[74|. The eQP approximation 
represents the diagonal propagator as a sum of two terms: 

= 4 - 5 (4n|4 + S A Il{S A ) . (75) 

The first describes decay processes, whereas the second 
can be associated with scattering processes. Inserting 
Eq. ([75]) into Eq. ([74")) we get a solution for the resummed 
Majorana propagator consistent up to the fourth order 
in the couplings: 

^ = e R [4 

- S R fl'^S A - \ {S R %S R + 5 A H|5 A )] Q A . (76) 
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The first term in the above formula describes Majorana 
decay, see Scc. lVI Cl whereas the remaining three terms 
describe the two-body scattering processes mediated by 
the Majorana neutrino. These are discussed in Sec. I VI D] 
Using definition of the retarded and advanced two- 
point functions, Eq. (|49p . and the Schwinger-Dyson equa- 
tion for the diagonal propagators, Eq. (|68|) . we find that 
the causal propagators in Eq. ()76[) are given by 



M-n 



R(A)) 



(77) 



Splitting the retarded and advanced self-energies into the 
vector and scalar components we can write the solution 
of Eq. ([771) in the form: 



S 



* 



R(A)) 



R(A) 



R(A)J 



(9-nV)) 2 -( M + n W 2 



ft h ± ? n p ' 



(78) 



where we have omitted flavor indices to shorten the no- 
tation and introduced 



n h =( q -ni) 2 -(M 



n 



From Eq. ([75)1 we can extract the spectral and hermitian 
propagators. To leading order in the Yukawas they read 



2 IX 



ft,, 



(79a) 
(79b) 



The on-shell condition is defined by ft h = 0. Expanding 
ft h to linear order in the Yukawas we find: 

fl h « {q 2 - A/ 2 ) (1 - 2r/7rM) 



v.jned 



+ (T/TrM)q 2 \n(\q 2 \/M 2 )-2qTl v h 

where jj^ ,med j s the medium-induced component of the 
hermitian self-energy in the on-shell rcnormalization 
scheme. In vacuum the on-shell condition is fulfilled 



for q 2 



M 2 , i.e. M is the physical vacuum mass. 



At non-zero temperatures the mass receives medium- 
induced corrections. To linear order in the Yukawas the 
effective mass is given by M 2 ss M 2 + 2qU^' med . For a 
hierarchical mass spectrum, which we consider here, the 
contributions of the hermitian self-energy are always neg- 
ligible and we will use ft h « q 2 — M 2 and fl h « d — M 
in the following. From Eq. (|T9"|) we can also deduce the 
effective width. To leading order in the Yukawas it is 
given by MT rj — qTl p . The minus sign in this definition 
ensures that the effective decay width is positive. One- 
loop contribution to the Majorana self-energy is derived 



in Appcndix[DJ In a CP-symmetric medium it is given 
by 



9w 
16tt 



[(A + ft)< J -Pi + (A t ft): i Pji]ip. (80) 



Therefore we can write the effective decay width in the 
form Tj = ■ (qLp/M 2 ), where Ti is the total vacuum 
decay width. For positive q 2 and q° the loop integral L p 
takes the form: 

L£ = 16tt J dnf : dU e p (2n) i 5(q-k-p)p tl 

x [1 + U(E k ) - f t {E v )\ . (81) 

For massless final states 2(qp) = M 2 . Therefore the defi- 
nition of the effective decay width inferred in Scc. lIIII from 
the requirement of successful MS subtraction is consis- 
tent with that implied by Eq. (j?9"jl . 

For the eQP Wightman propagators we can use the Ka- 
danoff-Baym ansatz. As can be inferred from Eq. (|75p. 
the corresponding spectral function reads 

J^-i^n^n^n^. (82) 

Substituting Eq. d78j into Eq. ((82]) wc obtain 



s p = - \ (n h tln h n 2 p + 4 n h #X n ? n h 



t v p t v p t v p nl)/{nl 



p 1 p 



(83) 



where we have again omitted the flavor indices. The 
second and the third terms in Eq. (|83|) vanish on the mass 
shell and can be neglected. Commuting fl h and Jt p in 
the first term and again neglecting contributions which 
are tiny on the mass shell we finally obtain for the eQP 
spectral function: 



[ft 



(84) 



Note that structures of Eqs. (JHJ) and (|79a[) are very simi- 
lar. Furthermore, as follows from Eq. (|11[) . in the limit of 
vanishing decay width both of them approach the delta- 
function. However, for a small but finite decay width the 
eQP spectral function is a better approximation to the 
delta- function than Eq. (|79ap . Therefore, we can approx- 
imate it by the usual expression, 

S p w (2tt) sign(q°)6(q 2 - M 2 )^ + M) , 

and at the same time keep finite-width terms in the di- 
agonal propagators. 



C. CP-violation in Majorana decay 

To go beyond the tree-level approximation and take 
into account CP- violating effects we need to consider con- 
tributions to the lepton self-energy that are of the fourth 
order in the Yukawa couplings. 
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One of them comes from expansion of the Majorana 
propagator in the one-loop self-energy. Substituting the 
decay term of Eq. (|76p into Eq. (|58[l we can write it in the 
form: 

Eg> (t, p) = - [ dU 4 dU 4 {2*) 4 5{q -k-p) (85) 

x (rfh) mn p R e%(t, q )s^(t,q)A § (t,k)eT(t,q)PL ■ 

Substituting Eq. (|8"5"j) and its CP conjugate into Eq. ([51]) 
we find that the resulting contribution to the divergence 
of the lepton-current has precisely the form (|66p . How- 
ever, the corresponding effective amplitudes are no longer 
equal: 



the masses. For a hierarchical neutrino mass spectrum 
the 'regulator' term is sub-dominant and this difference 
is numerically small. Note also that although Eq. ([87)) 
does not diverge in the limit of vanishing mass difference 
the approximations made in the course of its derivation 
are not applicable for a quasidegenerate mass spectrum 
2l| . For a consistent treatment of resonant enhancement 
within NEQFT we refer to [H|. 

The two-loop lepton self-energy is of the fourth order in 
the couplings to begin with. Therefore, for a hierarchical 
mass spectrum one can safely neglect the off-diagonal 
components of the Majorana propagators and replace 5? 
by the eQP one: 



xtr[0#(t, q){f + Mi)&t(t, q)P L pP R ] , (86a) 



'-'tdXr+Ni + ^(ch^Ni — 9w^2 mn (h h) m 



xtr[0#(t, q){i + M^Qtit, q)PLfP R } ■ (86b) 



.(2.1) 



(t,p) 



<m 4 dU 4 (2n) 4 6(p- 



q) 



x [{h^h)j 3 A 33 {t, q, k)P L CS^(t, q)P L A^(t, k) 

+ (h t h)%P R S%(t,q)CP R V jj (t,q,p)A^(t,k)] . (89) 



The matrices & R and A are evaluated on the mass shell 
of the i'th Majorana neutrino. The bar denotes CP-con- 
jugation and the trace is over Dirac indices. As compared 
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FIG. 7: Interference of tree- level and one-loop self-energy cor- 
rections. 

to tree-level result (frTT|) it additionally contains interfer- 
ence of the tree-level and one-loop self-energy contribu- 
tions to the Majorana decay amplitude, see Fig. [3 For a 
hierarchical mass spectrum we can use the approximation 

and a similar approximation for Q A . Using furthermore 
Eqs. ([79")) and (|80[) we find for the CP- violating parame- 
ter: 



AMfjMiTj 



{tfh)u{hMi) 33 (AM?-) 2 + (Tj/Mj qL f 
x^, 



(87) 



where p and q are on-shell momenta of the outgoing lep- 
ton and decaying Majorana neutrino respectively. In vac- 
uum = 0(q 2 ) sign(q )q )1 and the CP-violating param- 
eter takes the form: 



AMf-MiTj 



(tfh) u (tfh) 33 [AMI) 2 + (TjMf/MjY 



(88) 



The 'regulator' in the denominator of Eq. ([88)) differs 
from the result M,Tj found in [U by the ratio of 



Substituting Eq. ([59)1 and its CP conjugate into Eq. ([53]) 
we again find that the resulting contribution to the di- 
vergence of the lepton-current has the form (|66[) . The 
corresponding effective amplitudes read 



g w {tfh)% M t tr [A j3 (q, k)CP L fP R ] (90a) 
g w {tfh)% Mi tv[CV n (q, k)P L fP R \ , 
g w {h)h)% Mi tr [CV 33 (?, k)P L fP R ] (90b) 
g w {tfh)% Mi tx[A 33 {q, k)CP L fP R ] . 



They describe interference of the tree-level and one-loop 
vertex contributions to the Majorana decay amplitude, 
sec Fig.[5J To account for the contribution of the vertex 
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FIG. 8: Interference of tree-level and one-loop vertex correc- 
tions. 



correction to the decay width and the CP-violating pa- 
rameter we have to substitute the sum of and 
+ "^Ni^i^ an d a similar sum for the antiparti- 
cles into Eq. ([5Tj|) . The vertex contribution to the decay 
amplitude is of fourth order in the coupling and is neg- 
ligible compared to the tree-level term. Since we assume 
the medium to be almost CP-symmetric we can use, at 
leading order, CP-symmetric two-point functions in the 
loop integrals A 33 and V 33 . Then, at leading order in the 
Yukawa couplings, we find for the vertex contribution to 
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the CP-violating parameter: 
v _ ^3 {h)h% MiMj 



dK 



(91) 



^ (tfh)u qp 
x (2^) 4 % + ki + k 2 )(2ir) 4 5(k + k 2 - k 3 ){pk 2 ) 
x [A p (k 1 )S F (k2)S{'(k 3 ) + A F (fc 1 )S p (fc 2 )5^(fc 3 ) 
- A h (h)S p (k 2 )S j J(k 3 ) + A h (fc 1 )S F (fe)5^(fe 3 ) 
+ A^S^S^ih) + AMSMS^ih)] . 

The quasiparticle approximation and the KB-ansatz en- 
force two of the intermediate lines of the vertex loop to 
be on-shcll whereas the remaining line described by the 
hcrmitian part of the retarded and advanced propaga- 
tors remains off-shell. The three lines in square brack- 
ets in Eq. (f9"Tj) therefore correspond to different cuts 
through two of the three internal lines of the loop di- 
agram Fig. [I] (c) . Note also that only for one of the three 
internal lines the corresponding distribution function en- 
ters the result. 

The first possible cut described by the first line in 
square brackets corresponds to cutting the propagators 
of Higgs and lepton. One can interpret this cut as decay 
of the Majorana neutrino into a lcpton-Higgs pair which 
is followed by a subsequent t-channel scattering mediated 
by a virtual Majorana neutrino. Introducing 

K? (q, k) = 16vr / d< dU e p2 (2nf S(q -k 2 -p 2 )p^ 
x [1 + }^{E k2 ) - f t (E P2 )]M?St(k~ P2 ) , (92) 



we can rewrite the first term in Eq. (|9ip in a form which 
strongly resembles the form of the self-energy CP-viola- 
ting parameter: 



5 {h)h)l MjTj pK 3 
(Wh) vl {Wh) n Mf qp 



(93) 



In vacuum Kj can be computed explicitly and we recover 
the well-known result [j| : 



{h)h) 



[1 - (1 + Mf/Mf) In I 



M?/M$)]. (94) 

Adding up Eqs. ([55)) and (|9"4"|) we obtain the canonical ex- 
pression for the vacuum CP-violating parameter, Eq. (0. 

If the intermediate Majorana neutrino is much heavier 
than the decaying one then M^S\l « 1 and therefore 
Kj(q,k) pa L p (q). In this case we can also neglect the 
'regulator' term in the denominator of Eq. (|87|) . In this 
approximation the two CP- violating parameters have the 
same structure and their sum can be written in the form: 

, _ c vac PLp 
ti — tj 

qp 

Note that the combination of the distribution functions 
that enters the self-energy and vertex CP-violating pa- 
rameters, see Eqs. (|81[) and (|9"2"j) . is the same as that of 



ft<j> = 1 + Z^ 9 ~ ft 9 encountered in the derivation of the 
rate equations, see E q. ([321) . This result is in agreement 
with the findings of [STlIl El E3 usin S NEQFT and 
of [3l| based on imaginary-time thermal QFT. Note that 
older results featured a different dependence on the dis- 
tribution functions, with an additional term quadratic in 
the one-particle distribution functions which is absent in 
Eq. ([Ml) as well as in Eq. ([57]): 



1 + U ~ ft + 2 hh 



ft 



In [22| it was demonstrated that the result obtained us- 
ing thermal field theory can be reconciled with the re- 
sult of NEQFT calculation once causal Green's functions 
are used in the former. The two other cuts in Eq. ([9T]) 
are proportional to f N — f e and to f N + f, respectively. 
They vanish in the zero temperature limit and are usually 
Boltzmann-suppressed at finite temperatures, but can be 
relevant in specific cases [25| . 

The quantities that enter the rate equations are the de- 
cay, washout and CP-violating decay reaction densities. 
In the canonical approximation, i.e. when the quantum- 
statistical effects and effective masses of the Higgs and 
leptons are neglected, they are given by Eq. (|4"Tj) . If the 
thermal masses are neglected but the quantum-statistical 
effects are taken into account, there is an enhancement 
of the decay and washout reaction densities at high tem- 
perature, see Fig.[5J However, the inclusion of the ther- 




FIG. 9: Decay and CP-violating reaction densities with ther- 
mal lepton and Higgs masses, {X-y®.}, and with zero masses, 
{Xyjy.) m =o, for the two Majorana neutrinos Ni and N2. The 
values are normalized to the corresponding reaction density in 
the conventional approximation (Xj^i) ■ The thermal en- 
hancement due to quantum-statistical factors is overcompen- 
sated by the phase space suppression due to thermal masses 
at high temperatures. Note that we show only the self-energy 
contribution to the CP-violating reaction densities. 



mal masses turns this enhancement into a suppression at 
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high temperatures. It is explained by the decrease of the 
decay phase space. At intermediate temperatures the 
thermal masses become small relative to the Majorana 
mass and we observe a minor enhancement. For the CP- 
violating reaction density we observe a very similar be- 
havior. Given that for a hierarchical mass spectrum most 
of the asymmetry is typically generated by the lightest 
Majorana neutrino at zr ~ lnXi ~ O(l), where K\ is 
the washout parameter (see Appendix[E|, we expect the 
medium effects to induce a moderate enhancement of the 
total generated asymmetry. 



D. Majorana- mediated scattering 

Two-body scattering processes mediated by Majorana 
neutrinos violate lepton number by two units and play 
an important role in the washout of the generated asym- 
metry. In this section we derive the effective scattering 
amplitudes using NEQFT. This is an important part of 
our results. 



where we have defined: 



The last three terms in Eq. ([76)) contain the Wigner- 
transformcd one-loop Majorana self-energy: 

n| (t, g) = - 9w J dni dnfair^siq - P -k) 

x [(h''h) ij P L S^(t,p)P R A^(t,k) 

+ (h J 'h) ji P R PS^(t,p)PP L A^(t,k)] , (95) 



see Appendix lD 3 1 for more details. Combining them with 
Eq. (|58[) we find that their contribution to the divergence 
of the lepton-currcnt (|55p contains two Wightman prop- 
agators of leptons and two of the Higgs field. As we 
have argued above, these correspond to initial and fi- 
nal states in the kinetic equations. Therefore, we con- 
clude that these terms describe scattering processes de- 
picted in Fig. [2] and Fig. [3] As Higgs and leptons are 
maintained close to equilibrium we can safely use the 
Kadanoff-Baym ansatz for their propagators in the Ma- 
jorana self-energy Inserting Eq. (pj)) into the scat- 
tering terms of Eq. (|76|) we can then split the Majorana 
propagator into a lepton number conserving and lepton 
number violating part: 

^{l)=9 W J dUidU 4 p (27r) 4 8(q-p-k)A p (k)S p (p) 
x [(l-/f)(l + /|)^cfep) 

+ ( 96a ) 

y l i{q) = -9u, J dlli dU A p {2v)H(q -p- k)A p (k)S p (p) 

x [fIf^L j c(q,p) + f!f$^ v (q,p)] , (96b) 



lc - 
5 ij 



= (tfhiu [(1 - 6 ij )S R (q)P L fP R S 3 ^ (q) (97a) 
{Sl{q)P L fP R S»{q) +S i X{ q )P L pP R Sf{q))] , 



■ ' i 



LV - {h^h) Tl [(1 - 5^)S%{q)P R fP L S 3 i (q) (97b) 



+ ^(S%( q )P R pP L S^( q )+S2(q)P R pP L S^(q))] . 

Here we neglected higher order terms coming from the 
matrices Q R and 8a- The first terms in Eqs. (|97ap and 
(|97b[) corresponds to the second term in Eq. (|T6"|) . whereas 
the remaining terms correspond to the last two terms in 
Eq. ([To]) . Substituting the Majorana propagators (|9"6"1) 
into the lepton-currcnt (|55|) together with the one-loop 
lepton self-energy (|5"5|) we finally obtain 



J dn^ dn* x dn± 2 dn 



x A (3 (fc 1 )A p (fc 2 )S p (p 1 )S / ,(p 2 ) 



A LC {p 1 + k 1 , Pl ,p 2 )P l 



Plk 1 ;p 2 k2 



(98) 



Note that the zeroth component of the momenta in the 
above equation can have both signs. The effective ampli- 
tudes of the lepton number conserving and lepton number 
violating processes read 



LV 



LC 



(q,p 1 ,p 2 )=(h^h) Tl tT[^ v (q,p 2 )P L p i P R \ , (99a) 
(q,p 1 ,p 2 )=(h^h) Jl tl[y^ c (q,p 2 )P L f 1 P R ] . (99b) 



The functions A LC \q,pi,p 2 ) and A LV (q,pi,p 2 ) are sym- 
metric under the exchange of the momenta p\ and p 2 . 
This implies that the contribution of A LC in Eq. ([55)) 
vanishes. The terms of S^ LV diagonal in flavor space 
correspond to the RIS-propagator. Substituting Eq. ([75]) 
into Eq. (|99a[) and taking the trace we find that it con- 
tains only the scalar components of the retarded and ad- 
vanced propagators and is proportional to: 



2M 



2M' 



2 K-K 

(q 2 - M 2 ) 2 



{MTf 



[{q 2 - M 2 ) 2 + (A/r) 2 ] 2 



(100) 



Equation (jlOOp differs from the canonical result (TH} only 
in that the vacuum masses and decay widths are replaced 
by thermal ones. For a hierarchical mass spectrum this 
difference can be safely neglected. Introducing an ana- 
logue of the MS subtracted propagator, 

PiM = S l X(q)S%(q) - \^S p \q)S p \q) , (101) 

we can rewrite the lepton number violating effective am- 
plitude in a compact form: 

A LV (q, Pl , P2 ) = 2(rfh)l(p 1 p 2 )M i M j & ij (q) . 
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Next wc perform the trivial integrations over the fre- 
quencies using the Dirac-deltas in the quasiparticle spec- 
tral functions ([57)) and ([65]). Each Dirac-delta can be 
decomposed into two terms, one with positive and one 
with negative frequency. Therefore, the integration over 
the four frequencies gives rise to 2 4 terms, but only six 
of them satisfy energy conservation ensured by the re- 
maining delta-function. In a homogeneous and isotropic 
medium the one-particle distribution functions satisfy 

1 - fe(-p) = flip) , 1 + U(-k) = -/#) , (102) 

and the diagonal Majorana propagators have the prop- 
erties 

S i ;(-q) = -S i j(q), S%(-q) = S${q). (103) 

Upon substitution of the resulting self-energy into 
Eq. ([55) and the use of Eqs. ([TU^]) and ([TUB")) the remain- 
ing six contributions in the lepton-current can be conve- 
niently written as 



dll 



p\p 2 k^k 2 



■pkik 2 ;piP2 



— (txt) 



-pPlp 2 \kik 2 



p/ CT (sXs) , —(tXt) ■> 



-pPlk 1 ;p 2 k 2 



(104) 





sxs) 








(txi) 








(txt) 




tl<r¥4> 




(txt) 







where we have defined the effective scattering ampli- 
tudes: 

20* (piPa) ^irh r r \l,.\l r A J ui. ! , (105a) 

25* ipm) *l(h!h)? j M i M J & ij (q t ) , (105b) 

25* (P1P2) ^(/^^M^-fe) . (105d) 

The momenta of the Majorana neutrinos are related to 
the momenta of the initial and final states by q s = Pi + k\ 
and q t = pi — k 2 ■ From Eq. (|105[) we see that the obtained 
amplitudes contain only sxs and txt interference terms. 
Indeed, in the products of the Majorana propagators in 
Eq. (jlOip both of them depend on the same momentum. 

The missing cross terms emerge from the two-loop (ver- 
tex) contribution to the lepton self-energy. As we have 
mentioned in Sec.[Vl within the discussed assumptions 
and approximations the second and third terms on the 
right-hand side of Eq. ([55]) describe scattering processes. 
Since they are of the fourth order in the Yukawas, we 
can replace the full Majorana propagators by the diag- 
onal propagators. The third term, Eq. (|62[) . corresponds 
to lepton number conserving processes and does not need 
to be discussed further. The second one, Eq. (pxTj) , is given 
by 

£g' 2) (*>*>)= / d K dUi i dUi 2 (2n) i 5(p + k 1 -p 2 ~k 2 ) 
x {h)h)l [P R S%(t,p 2 + k 2 )CP R S^(t, - P2 )P L 
xCS i i(t,p 2 -k 1 )P L A^(t,k 1 )A^(t,-k 2 )] . (106) 



We substitute Eq. (|106p into the equation for the lepton- 
current (|55p and perform the steps preceding Eq. (|104p . 
Using furthermore relations f|102[) we find 



= - dn 



-iPlkip 2 k 2 



(sxt) -pp\kr,p 2 k 2 



2S 



_,(uxt) -pp 1 p 2 ;k 1 k 2 



r-,(iiXt) -pk 1 k 2 ;p 1 p 2 
"Weft W«« 



where we have introduced 

= 4g w (piP2) £ U.U/J? (ffh)l 

xU{S% t )S j i(q s )}, (107a) 
x^(rth)* ij S i X(st t )S%(q u )}, (107b) 

xmrthf^SKqJS^qt)}, (107c) 

and q u = p 2 — k 2 . The combinations of the momenta 
appearing in the products of the Majorana propagators 
clearly indicate that the above amplitudes correspond to 
the interference terms of the s-, t-, and u-channel contri- 
butions. 

Combining Eqs. (| 105[) and (|107f> we obtain for the ef- 
fective amplitude of ££ O 4>(f> scattering: 

E u ^ = 2(p lP2 )J2(tfh)lM i M j 
x [4^- (ft) 

+ Si(q t )S»(q u )+S i i(q u )S j i(q t )] , (108) 

whereas the effective amplitude of £<p <-> £<fi scattering 
reads 

= 4(pi P2 ) £ U(h^h)ijMiMj 
x [2^> ij (q 8 )+2^ ij (q t ) 

To compare the obtained expressions with the vacuum 
results of the canonical computation Eqs. ([T6|) and (jTTJ) , 
we evaluate the retarded and advanced propagators at 
zero temperature. In vacuum L p (t,q) = 0(q 2 ) sign(<7°) fl. 
Therefore for positive qo it follows from Eq. (|77[) that 



(?) 



[q 2 -M?±l6(q 2 )r t /M iq 2 



(109) 



In the vicinity of the mass shell of the respective Majo- 
rana neutrino q 2 = M 2 and we find 



SI = -Pi 



si = -p; 



(110) 



For t- and u-channels the imaginary parts of Eqs. (|10p 
and (|109[) vanish, so that in the vacuum limit 

&> ij (q)=P*(q 2 )P j (q 2 ), 
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and is symmetric with respect to j o j. The 

squares of the t- and u-channel propagators in Eq. (TIT]) 
give identical contributions to the reduced cross section 
of ££ O- (j>(f> process. Therefore, upon substitution of 
Eq. (|108p into Eq. ([45]) and the use of the i «-> j symme- 
try, we recover for the reduced cross section the canonical 
result ([44]). Relations (|110|) imply that the off-diagonal 
components of Pjj and coincide. For the s-channcl 
the diagonal components of are given in the vacuum 
limit by 

^ n{qs> [(a-Mfy + pi/Mis)*]*' 

Thus, in the vicinity of the mass shell, s sa Mf, the 
canonical expression for the RIS-subtracted propagator, 
Eq. (|14[) . coincides with the expression obtained from first 
principles. This justifies results of the earlier calcula- 
tions. For the reduced cross section of t<j> £(f> scattering 
we recover Eq. ([4"5l) . 

At finite temperatures L p (t,<ft) is not zero even for t- 
and it-channels. In other words, the medium effects in- 
duce additional contributions to the effective decay am- 
plitudes. However, these contributions are proportional 
to the coefficients c,-. Numerical analysis shows that for 
the two chosen sets of parameters, see Appcndix[El the 
additional correction typically do not have any sizable 
impact on the reaction densities. 

A quantity relevant for the numerical analysis is the 
ratio z{"ff!?) /Hs. The dependence of this ratio on the 
dimensionless inverse temperature is presented in Fig. 1101 
and Fig. [11] If the approximate expression ([42]) is used, 
then the reaction density of £<j> £(j) scattering becomes 
negative for 2 < z < 3 for the first set of the param- 
eters whereas for the second set of the parameters it 
turns negative for 0.5 < z < 1. A qualitatively simi- 
lar behavior has also been observed in [4l|. This rather 
counter-intuitive result can be traced back to the behav- 
ior of the RIS part of the effective amplitude (|16[) which 
is negative in the vicinity of the mass shell. Its sign is 
not fixed by physical requirements, since it constitutes a 
sub-leading contribution to the washout rate. As can be 
inferred from Fig. 1101 for the first set of model parame- 
ters the quantum-statistical corrections render the RIS 
subtracted reaction density of £(f> O 1(f) scattering posi- 
tive in the whole range of temperatures. However, this 
is merely a numerical coincidence. For the second set of 
parameters, see Fig. [11] the reaction density of £(j> <-> £<fi 
scattering remains negative for 0.4 < z < 0.9. It is nev- 
ertheless important to note that even in the region where 
the reaction density is negative the quantum-statistical 
corrections shift it upwards as compared to the result of 
the canonical computation. 

As far as ££ -H> <f>4> scattering is concerned, the 
quantum-statistical corrections enhance the reduced re- 
action density at high temperatures by about 50%. As 
the temperature decreases, the reaction density com- 
puted using Eq. (|35[) approaches the one computed using 
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FIG. 10: Washout reaction densities due to ££ o <fxj> and 
lxj> f-> £<f> scattering processes for benchmark point 1 in the 
approximation of massless leptons and Higgs. Shown are the 
reaction densities computed using the Boltzmann approxima- 
tion (I42|) (thin lines), and taking into account the quantum- 
statistical effects, (|35|) (thick lines). The (RIS subtracted) 
reaction densities may be negative as they are not the 

physical rates for the 2ft2 scattering process. 
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FIG. 11: Washout reaction densities due to ££ o <f>4> and 
£tj) o £<j> scattering processes for benchmark point 2. Com- 
pare Fig. [10] 

Eq. (|4"2"]) as one would expect. A similar behavior is also 
observed for the reaction density of £<fs <H- £<p scattering. 

In Figs.[TTJl and [TT1 we neglected the thermal masses of 
initial and final states. To estimate the size of the mass 
corrections, in Fig.[T2] we plot ratio of the reaction den- 
sity of the ££ -f-> <fxj) scattering computed using thermal 
masses of the Higgs and leptons with and without the 
quantum-statistical terms to the canonical one for the 
two sets of parameters. If the quantum-statistical effects 
are neglected, the thermal masses lead to a ~ 15% sup- 
pression of the reaction densities. On the other hand, an 
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FIG. 12: Ratio of the reaction density computed using ther- 
mal masses with (thick lines) and without (thin lines) the 
quantum-statistical terms to the canonical one for the two 
sets of parameters. 



enhancement induced by the quantum-statistical terms 
to a large extent compensates the mass-induced suppres- 
sion. As a result, the deviation from the canonical reac- 
tion density does not exceed ~ 5% in the whole range of 
temperatures. 

To compare the relative importance of the (inverse) 
decay and scattering processes in Fig.[T3] and Fig.[H] we 
also present the uniformly normalized reaction densities. 
In both cases we observe a qualitatively similar picture: 

<7>/(2riMj), BM1 




FIG. 13: Washout reaction densities due to it <-> <f)<f> and 
£<f> f-> l(j> scattering processes for benchmark point 1 in the 
approximation of massless leptons and Higgs. For comparison 
the washout reaction densities for Ni and N2 (inverse) decays 
are shown as well. Note that the normalization differs from 
the one used in Fig. [10] and Fig. 1111 For the present choice of 
parameters contributions by 2 <->• 2 scatterings are strongly 
suppressed by the smallness of the couplings. 



for the chosen sets of parameters reaction densities of 
the scattering processes are strongly suppressed by the 
smallness of the Yukawa couplings as compared to those 
of the decay processes. 

{j)/(2T 1 Mf), BM2 
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FIG. 14: Washout reaction densities due to It o (j><f> and 
£(j> <H> £(j> scattering processes as well as (inverse) decays for 
benchmark point 2. Compare Fig. 1131 

In the rate equation for the lepton asymmetry, see 
Eq. (|29a|) . the total washout rate is given by a sum of 
reaction densities for decays and scattering. Whereas 
(7]^) and (7^) are positive, the RIS subtracted reac- 
tion density (7^ ) can be negative at some temperatures. 
If the total washout rate would turn negative, it would 
lead to a spurious self-enhanced generation of the asym- 
metry. In Figs.[T3] and [T?] we show the sum of the re- 
actions densities. For both parameter sets it is positive. 
This sum should always be positive. If the quantum- 
statistical terms are neglected this can be demonstrated 
explicitly. In the absence of the quantum-statistical cor- 
rections, the results of this section revert to the ones 
discussed in SecHU Using the expression for the RIS- 
propagator, Eq. (fl~5|) . we can rewrite the RIS subtracted 
scattering amplitude (|16p as a difference of the unsub- 
tracted one and the RIS term. Since the latter is propor- 
tional to 6(s — Mf) the integration in Eqs. ([4"3")l and 
is trivial and we obtain after some algebra: 



an) 



The total washout rate is then given by 4(7~) + 4(7~) 
and is positive as a sum of two positive functions. In 
other words, even though (7ji ) can be negative at some 
temperatures, the total washout rate is always positive. 
The form of the unsubtracted scattering reaction density 
that can be inferred from Eq. (|111[) is another manifesta- 
tion of the double-counting. If we had not subtracted the 
RIS contribution, we would have counted contributions 
of the inverse decay processes twice and ended up with 
an incorrect prediction for the generated asymmetry. 
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VII. HIGGS CONTRIBUTION 

In the preceding section we have approximately taken 
the gauge interactions into account in the form of ef- 
fective masses of the Higgs and leptons. 10 The ther- 
mal masses arc of order of gT and large enough to 
influence the values of the reaction densities quantita- 
tively. In particular, the Majorana neutrino decay can 
become kinematically forbidden when the sum of the 
masses of lepton and Higgs exceeds the heavy Majo- 
rana neutrino mass. For even higher temperatures (with 
m^(T) > mi(T) + Mi(T)) the Higgs decay channel into 
a lepton-Majorana pair becomes kinematically allowed 
instead and can contribute to the asymmetry since it vi- 
olates CP. For simplicity we do not take modified disper- 
sion relations into account here, see [HI HH, HH, H3, E3 > 
but use the simple picture of temperature dependent 
thermal masses as an estimate: 



—92 + —9 n 
16 y2 16 y 

3 2 . 1/2 

32 52 + 32 5 



4 * 



2 



(112) 



where we use the temperature dependent values of the 
SU(2)l, U(l), top Yukawa and Higgs self-coup ling s 52, 
g' , h\ and A assuming a Higgs mass of 115 GeV [76|]. We 
also ignore that the thermal mass of leptons might be 
better approximated by the 'asymptotic thermal mass' 
\2m(. in a kinematic regime in which their momentum 
is such that p 2 ~ (gT) 2 [Z^. We do also not take into 
account, in our quantitative analysis, the thermal cor- 
rections to the Majorana neutrino masses as they are 
negligible compared to their vacuum masses. 11 





P 



FIG. 15: Tree level contribution and one-loop corrections to 
the (anti-)Higgs decay amplitude. The additional arrows il- 
lustrate the direction of momentum flow. 



10 We will not attempt a fully consistent inclusion of gauge inter- 
actions here. Their effect has been addressed systematically for 
the production rate of Majorana neutrinos in various tempera- 
ture regimes in [34, 35, 37, 7jJ, but not for CP- violating rates up 
to now. In this work we model gauge interactions in a simplified 
way by including thermal masses, similar as has been done pre- 
viously in the context of thermal field theory computations [id . 
In this way the results obtained within NEQFT can be compared 
to previous computations, and may ultimately be compared to a 
full treatment of gauge interactions within NEQFT. 

11 Note that they may be relevant in the different context of reso- 
nant leptogenesis if the size of the thermal correction is compa- 
rable to the mass splittings | Mi — Mj | as they can influence the 
resonance in this case. 



The CP- violating decay of the Higgs arises from the in- 
terference of the tree-level, self-energy and vertex graphs 
depicted in Fig. [15] Although in this case the decaying 
particle - the Higgs doublet - is very close to thermal 
equilibrium due to the Yukawa and gauge interactions 
of the Standard Model, the Majorana neutrino in the fi- 
nal state may deviate from equilibrium, so that the third 
Sakharov condition is fulfilled. 

Using the expression for the divergence of the lepton- 
current, Eq. (|55[) . we can extract the corresponding CP- 
violating parameter. To calculate the self-energy con- 
tribution it is convenient to rewrite the one-loop lepton 
self-energy (|5"5]) in the form: 

x &^(t,k)P R C{y*) T {t,q)C- l P L , (113) 

where the transposition is only in Dirac space and we 
have used one of the properties of Majorana propagator: 

(^)(X,q) = C{^) T {X,-q)C- 1 . (114) 

In its CP conjugate the Yukawa couplings are replaced 
by their complex conjugates and the propagators by the 
CP conjugate ones. 

Similar to the case of the Majorana decay, substituting 
Eq. (|113p into Eq. ([55]) we can define effective Higgs decay 
amplitudes: 



. = 9w y^(fe t fe)r 



mn 

x tr + Mi)e^(q)P R fP L ] , (115a) 

inn 

mn 

x tr [eS*(?)(0 + Mi)e%(q)P R pP L ] . (115b) 



The overall factor g w in Eqs. (|115[) comes from summa- 
tion over the doublet components of the (decaying) Higgs 
particle. To leading order in the couplings: 

3 ^J\r, ~ ^9w[(h f h), H (p-q) 

+ faQ(tfh)$ j M i M j S»(pL p )] , (116a) 
Vw, ~ 2 9w [(h f h) u (p ■ q) 

- fe$i(rthfi j M i M j S i h i (pL l ,)] , (116b) 

which, up to the relative sign in square brackets, coin- 
cides with the amplitudes (156"]) . The corresponding CP- 
violating parameter reads 



is 



MiTj pL p 



(h^h)U^h) n M 2 - M 2 pq ' 



(117) 



where p and q are the momenta of the on-shell final lep- 
ton and Majorana neutrino with positive zeroth compo- 
nents respectively. The direction of momentum flow is 
as defined in Fig.flol Although L p in Eqs. ([57]) and (flTT[) 
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is one and the same function, because of the different 
kinematic regimes the explicit result in terms of the dis- 
tribution functions differs for the Higgs decay: 



L p (t,q) = 16ir J dnfdllfair^Siq+p 



x[U(E k ) + f e (E p )} 



(118) 



see Appcndix lD 31 Note that our result is different 
from the one presented in [l(], El, [z3- Instead of the 
fi — f e — 2f^f e dependence, it is proportional to a sum, 
f<t> + fti °f the two distribution functions. This depen- 
dence can also be obtained in the framework of real time 
thermal field theory using causal n-point functions, com- 
pare [H| . The derivation within the Kadanoff-Baym for- 
malism gives certainty concerning the sign of the con- 
tribution by Higgs decay. The CP-violating parameter 
(|117p has an opposite sign relative to that for Majorana 
neutrino decay. However, it is canceled by the relative 
sign in Eq. (|126p . 

To calculate the vertex contribution we use Eq. (|114|) 
and represent the two-loop self-energy (1891 ) in the form 

4 2 ' 1) (*>f) = / dl^dUt (2nf6(p + q - k) (119) 
x [{h)h)%A n (t, -q, ~k)P L (S^) T (t,q)CP L A^(t, k) 
+ (h^h)%P R C(S^) T (t, q)P R V jj (t, -q, -p)A^(t, k)] . 

Its CP conjugate again differs by the conjugation of the 
Yukawas and propagators. Substituting Eq. (|119p and its 
CP conjugate into Eq. ([5"5"j) we obtain for the correspond- 
ing effective amplitudes: 

^ M = 9w {h)h% M^A^-q-^CPLfPn] 
+ g w {h)h)% Mi tr[C%(- 9 , -k)P L fP R ] , (120a) 

^lN t = 9^h)l M^CV^-q-^PLfPn] 
+ g w (h f h)% M<tr[A i3 -(-g, -k)CP L fP R ] . (120b) 

Just like for the self-energy contribution we observe that 
the overall sign of the vertex contribution to the Higgs 
decay amplitude is opposite to that in the Majorana de- 
cay, compare Eas. (|90|) and (|120|) . The corresponding 
CP- violating parameter reads 

x {2n) 4 S{q - hi - k 2 ){2ir) 4 S(k - fc 2 - k 3 )(pk 2 ) 
x [A p (h)S F (k 2 )S{'(k 3 ) + A F (fc 1 )S p (fc 2 )5i J (fc 3 ) 
- A,(fc 1 )S p (fc 2 )5^(fc 3 ) - A h (h)S F (k 2 )S^(k 3 ) 



vertex CP-violating parameter depends on the tempera- 
ture as well as on masses of the Majorana neutrinos. For 
definiteness, let us assume a strongly hierarchical mass 
spectrum, Mj 3> > Mi. In this case contribution 
of the last two lines in Eq. (|121[) is strongly suppressed. 
Integrating out the delta-functions we find for the con- 
tribution of the first cut: 



v 



lJ2 



2 ^ {Wh) u {h^h) 3J M? qp ' 
where the loop function Kj is now defined as: 



(122) 



K?(q, k) = 16tt J dUl dH e p2 (27r) 4 5(g + p 2 - k 2 )p^ 2 

x [U(E k2 ) + f e (E P2 )]M?St(k-p 2 ). (123) 

Note that the vertex CP-violating parameter (|122[) has 
an opposite sign relative to that for Majorana neutrino 
decay. Similarly to the self-energy contribution we ob- 
serve that the 1 — + combination is replaced in 
(|123p by f e + f^. For a milder mass hierarchy the two 
other cuts can become important. Their contributions 
are proportional to 1 — f t — f N and + f N respectively. 

The first-principle computation gives for the Higgs de- 
cay contribution to the evolution of the lepton current 
an expression similar to Eq. (|24|) : 



dYL 



dz 



z 



dn 



qkp 
-qp;k 



-pqp;k 



(124) 



We do not discuss 'extra' terms here which would arise 
from 'naive' Boltzmann equations. To write this as a 
rate equation we need to repeat the steps in Sec. lIVI 
which lead to Eq. (f3"6")) . For a general process aN <R- b 
(where we allow for deviations from equilibrium in f N ) 
we have Eq. (|C3[) . Therefore we obtain, for Higgs de- 
cay, the following contributions to the rate equations for 
lepton number and Majorana neutrino abundance (see 
AppendixIC]): 



sU dY Lt 



z dz 
sH dY Ll 



£d>A „vq 70, i I 



f. 



(125a) 



z dz 

dYL, 



H I /-i req req*. ^/w, D \ 



sH 



dz 



W,, . x/( 1_ /iV») w 



(125b) 
(125c) 



A p (ki)S h (k 2 )Sp (k 3 ) — A F (ki)S h (k 2 )S p J (fc 3 )] . We introduced the decay and washout reaction densities, 



The last three lines of Eq. (|12ip correspond to the three 
possible cuts of the vertex graph. Similarly to the Ma- 
jorana decay only two of the intermediate states can be 
on-shell and for only one of them the corresponding dis- 
tribution function enters the result. The value of the 



(XrfJ and {Xy% t }, for Higgs decay: 



(X^)=J dUZ k i<p (27rrS(q+p-k)XE^J l4> , 
(X^={X(l-f%)tfJ, (126) 
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where now = (f^ 9 + f^ 9 ). In complete analogy to 
Eq. (|30p the total amplitude and the CP-violating pa- 
rameter for (anti-)Higgs decay are defined as 



e <t>,i 



2S 



Similarly to the Majorana neutrino decays we also define 
an averaged CP-violating parameter as 



<7 D 



By comparing Eqs. (|125a|) and (j36)) we observe that, ig- 
noring Af N .fig/T contributions, the difference to the Ma- 
jorana neutrino decay contributions amounts to the re- 
placements (X 7 g) -> -{XrfJ and (A 7 J£) -> (X 7 ^). 
We therefore obtain Eq. (|29ap with an opposite sign for 
the CP-violating source term. This sign cancels the rel- 
ative sign of the CP-violating parameter such that Ma- 
jorana neutrino decay and Higgs decay contribute effec- 
tively with same sign. Similarly, for the contribution to 
the Majorana neutrino rate equation: 



sH dY Ll 



z dz 
sH dY Ll 



z dz 
dYh 
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sH 
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req.^hk D 
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To compute the thermally averaged CP-violating param- 
eter we take into account the temperature dependent evo- 
lution of the lepton and Higgs masses (|112p jlQj . The av- 
eraged CP- violating parameter in the Higgs decay and in 
the Majorana decay as functions of the inverse temper- 
ature are presented in Fig. [TBI and Fig.JTT] At very high 
temperatures the magnitude of (ej) can be much larger 
for Higgs decay. As the temperature decreases the CP- 
violating parameter for the Higgs decay approaches zero. 
This is explained by the shrinking of the available phase 
space in the loop integrals (|118[) and (|123p . 
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FIG. 16: Averaged self-energy CP-violating parameters for 
Majorana neutrino and Higgs decay for benchmark point 1 as 
a function of the inverse temperature. Thin lines represent 
the value in the zero temperature limit. With conventional 
dispersion relations, the decay Ni — > £(j> is active at z > 4.93 • 
10" 1 but replaced at high temperatures (z < 1.18 • 10" 1 ) by 

but 



" but replaced at high temperatures (z < 1.18 • 10 
Nit The decay N 2 -> l(j> is active at z > 1.54 ■ 10" 



replaced at high temperatures (z < 3.62 -10 ) by <f> — > A^. 
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FIG. 17: Averaged self-energy CP-violating parameters for 
Majorana neutrino and Higgs decay for benchmark point 2. 
See caption of Fig. 1161 



VIII. CONCLUSIONS 

In this paper we have studied leptogenesis in the type-I 
seesaw extension of the Standard Model using the 2PI- 
formalism of non-equilibrium quantum field theory. 

The asymmetry generation can, in the case of thermal 
leptogenesis, be approximately described by rate equa- 
tions. Usually these statistical equations are treated as 
a 'black-box' in the sense that their form is assumed 
given and model specific amplitudes are inserted by hand. 
Indeed, this approach is supported by the observation 



that it describes the free decay in the zero temperature 
limit correctly and inherent inconsistencies (namely the 
'double-counting problem') can be resolved in exact equi- 
librium. However, out of equilibrium it is not obvious 
whether the subtraction of real intermediate states works 
to all orders and how amplitudes computed in thermal 
field theory enter kinetic equations. These issues can 
be completely avoided in a systematic treatment within 
non-equilibrium quantum field theory. 

Only in recent years fundamental questions related to 
the non-equilibrium statistical description received more 
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attention. The progress here is mainly based on the 
2PI-formalism which is known to yield consistent quan- 
tum kinetic equations without double-counting. These 
equations can be reduced to a system of Boltzmann- 
like kinetic equations for quasiparticles which can easily 
be compared to the conventional results. In the course 
of the derivation necessary approximations and the re- 
lated physical assumptions have to be specified explicitly. 
Therefore, this approach enables a deeper insight into the 
dynamics of the asymmetry generation. 

In the conventional analysis the minimal set of in- 
teractions is obtained at order 0(h 4 ) of the perturba- 
tive expansion. We have complemented existing analyses 
based on the 2PI-formalism by the computation of fur- 
ther processes which appear at this order. Starting from 
a system of Kadanoff-Baym and (equivalent) Schwinger- 
Dyson equations for leptons and heavy Majorana neu- 
trinos we have derived Boltzmann-like quantum-kinetic 
equations for the lepton asymmetry. They include (in- 
verse) decays of the heavy Majorana neutrinos as well 
as two-body scattering processes mediated by the heavy 
neutrinos: 
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Because all terms in this equation are proportional to J-, 
a combination of the distribution functions which van- 
ishes in equilibrium, the obtained equations are free of 
the double-counting problem and no need for the real in- 
termediate state subtraction arises. Together with the 
systematic derivation of the effective decay and scatter- 
ing amplitudes S this is the main result of the present 
work. The individual amplitudes arise as combinations 
of different 2PI contributions. The vertex contribution 
to the CP-violating decay amplitude is obtained as cut 
of the 2PI 'mercedes' diagram. The same graph yields 
also the s x t contributions to £<p <->■ £<fi and the uxt con- 
tribution to ££ -H> 4>(f> process. To extract the self-energy 
contribution, the off-diagonal elements of the Majorana 
neutrino propagator have to be taken into account. In 
addition, an extended quasiparticlc approximation needs 
to be employed in order to obtain the s x s and t x t 
contributions to £(p £(f> as well asuxii and t x t con- 
tributions to ££ o 4>(j> scattering from the 'setting-sun' 
diagram. 

In the zero temperature limit the effective amplitudes 
reduce to the canonical ones. In particular, the form 
of the resulting amplitudes for £<fi <-> £<fi scattering co- 
incides with the RIS subtracted amplitudes encountered 



in existing calculations. At finite temperatures the ef- 
fective amplitudes receive thermal corrections. Medium 
corrections to the Majorana decay amplitudes into lep- 
tons and antileptons are 0(h 4 ). They are small compared 
to the tree-level vacuum contribution and are therefore 
negligible for the total decay width. On the other hand, 
they play an important role for the CP- violating source- 
terms, which are proportional to the difference of the 
two amplitudes. We find that medium corrections to the 
CP-violating parameter are linear in the particle num- 
ber densities. Although there is a partial cancellation of 
the bosonic and fermionic contribution, the CP-violating 
parameter is enhanced. In the effective scattering ampli- 
tudes the medium corrections affect only the 'regulator' 
term in the denominator of the Breit-Wigner propaga- 
tors. Due to the smallness of the Majorana decay width, 
which is constrained by the light neutrino masses, nu- 
merically these corrections are very small in the case of 
non-degenerate Majorana neutrinos. 

Taking SM interactions into account in the form of 
thermal lepton and Higgs masses results in a suppression 
of the phase space for the Majorana neutrino decay and 
the enhancement of the CP- violating parameters is over- 
compensated. At even higher temperatures, when the 
effective Higgs mass exceeds the Majorana masses, the 
CP-violating decay of the Higgs into a lepton-Majorana 
pair can become kincmatically allowed instead. At these 
temperatures the averaged CP-violating parameters for 
Higgs decay exceeds that obtained for Majorana decay in 
vacuum by orders of magnitude. The signs of the corre- 
sponding CP- violating parameters are opposite but their 
contribution to the lepton asymmetry has the same sign 
(at least in the limit of hierarchical Majorana masses). 
These results are in qualitative agreement with earlier 
studies based on thermal field theory and may ultimately 
be compared to a full treatment of CP- violating decays 
out of equilibrium, including gauge interactions. 

We have also derived the corresponding rate equations 
for abundances of the participating species. They are 
obtained as expansion in small deviations from equilib- 
rium (fi/T and Af N .) and represent the hydrodynam- 
ical approximation of the Boltzmann kinetic equations. 
As compared to the standard (zero temperature) result 
they are improved in that the obtained coefficients in- 
clude medium corrections to the quasiparticle properties 
and take into account quantum-statistical effects. We 
compare with the result obtained if the amplitudes are 
computed in thermal quantum field theory and the RIS 
subtraction is performed manually. We find that there 
are differences at higher order in the expansion param- 
eters. The coefficients - reaction densities - reflect the 
interplay between the medium enhancement of the effec- 
tive amplitudes and the phase space suppression induced 
by the thermal masses of Higgs and leptons. At very 
low temperatures the reaction densities approach their 
canonical limit. 

Since for a hierarchical mass spectrum most of the 
asymmetry is typically generated by the lightest Majo- 
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rana neutrino at temperatures of the order and smaller 
than its mass, we expect a moderate enhancement of the 
total generated asymmetry is possible for a typical av- 
erage to strong washout scenario. For a detailed phe- 
nomenological analysis it is necessary to include further 
phenomena such as flavour effects and AL = I scattering 
processes which contribute to the washout at 0(g 2 h 2 ). 
Additional quantum effects beyond the present analysis 
are relevant for non-standard scenarios in which the Ma- 
jorana neutrinos have degenerate masses or if they are 
not as close to thermal equilibrium. 
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Appendix A: Kinematics 

The reaction densities contain distribution functions 
of the initial and final states, which depend on the cor- 
responding energies. Therefore to compute the reaction 
densities we need to analyze the kinematics of the decay 
and scattering processes. 



1. Decay 

To compute the CP- violating reaction density and the 
decay reaction density we need to evaluate the integral: 

(X 7 ^) = J dnf/pj dILtdILi(2n)*6(k+p-q) 

xXE N [l-f? + f?]. 

For the washout reaction densities or reaction densities 
for Higgs decay we have similar expressions. The inte- 
gration over d 3 k can be performed trivially and yields 
k = q — p. Using |fc| = (q 2 + p 2 — \q\\p\ cosB)* and inte- 
grating over 9 we remove the remaining Dirac-dclta and 
obtain 

E + 2ir 

J dUidU^rSik+p-q) ^^JdE p J d £. 

E- 

(Al) 

The integration limits are given by 

E± = § [E q (1 + xi - x ) ± \q\\? (1, xi, X<t> )] , (A2) 



where xe = m 2 /q 2 , x^ = mhjq 2 and A(x, y, z) = x 2 +y 2 + 
z 2 — 2xy — 2xz — 2yz is the usual kinematical function. 
For an on-shell heavy neutrino q 2 = M 2 . If me = = 
then A(l,x^,X0) = 1 and the above expression simplifies 
to E± = \{E q ± \q\). On the other hand, if M = mi + 
?7i0 then \(l,xe,x$) = and therefore E + = E~ = 
E q (me/M). Since the integration limits coincide in this 
case, the integral vanishes. 

Combining it with the integration over d 3 q and using 
the isotropy of the medium we find 

oo -E+ 

<*7&) = ^hJ dE i fN J dE p XE N [l f? + . 

M E- 



If quantum-statistical effects are neglected then both the 
CP-violating parameter and the total tree-level decay 
amplitude are momentum independent and the integra- 
tion can be performed analytically. We get: 

( X « ^XMfT.TK, (f) , 

where g^ = 2 is the number of the Majorana spin degrees 
of freedom. Similar results are obtained for Higgs decay 
and washout reaction densities. 



2. Two-body scattering 

For 2f>2 scattering processes the reaction density is 
defined by 

( 7 #) = J dSf$"» (2n)*5( Pa + Pb - V i - Vi) 

x f?f?(l ± /f )(1 ± /f )E ab ^ ij . (A3) 

To reduce it to a form suitable for the numerical analysis 
we insert an identity: 

1 = J ds J d i qS(p a + p b - q)5 + (q 2 - s) , 

into Eq. (|A3|) . The resulting expression can be inter- 
preted as a product of the inverse decay and decay am- 
plitudes integrated over the 'mass' and energy of the in- 
termediate state: 

-+ 1 ds J 'w? 5+(q2 s) 

x J dn; a dn b Pb (2ir) 4 6(p a + Pb -q) 

x / dW pt dW {2^f5{q-p t - P] ).{M) 



29 



For the third line we will use Eq. (|A1[) . For the last line it 
is more convenient to us a different representation. Inte- 
grating out the delta-function we obtain for the last line 
in Eq. ([M]) : 



PI 



in J 4tt q a \pi\ - Ei\q\ cos8 



(A5) 



Note that not all angles are kinematically allowed, see 
Eq. (|A8p below. As a product of Lorentz-invariant quan- 
tities the integral is also Lorentz-invariant. We can there- 
fore boost to the center-of-mass frame where q = and 
Qo = \/s- By energy-momentum conservation \pi\/qo = 
h\z (1, Xi, Xj), where q 2 = s now. The angle integra- 
tion can be partially reduced to integration over the 
Mandclstam variable t = (p a — pi) 2 using the relation 
dt = 2|p a ||pi|c?cos6' a i. Note that the azimuthal angle ip a i 
is Lorentz-invariant by itself. Therefore, boosting back to 
the rest frame of the medium we can write the left-hand 
side of Eq. (|A5|) in the form 



\ 1 f dl Pai 

A 2{l,x a ,x b ) ■ — / — — 
»7r / 2ir 



dt 
s 



Integrating furthermore over dq° and using the fact the 
the integrand is independent of the orientation of q, we 
finally obtain 



oo oo 



(7. 



ab 



l 



(A6) 
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x\-Hl,x a ,x b ) I dE a f?f? 



lit V 



o t- 

where \fs min = max(X] m. lmt , J2 m fin)- Just like for 
particle decay, the integration limits E^ are given by 
Eq. (|A2|) but with xg and x^ replaced by x a and x b respec- 
tively and the three-momentum given by \q\ = (E 2 — s) = . 
The range of integration over t is given by 



- -[(1 + X a - X b )(l +Xt- Xj) 

=F \?(l,x a ,x b )\i(l,Xi,Xj)] . 



In particular, for massless initial and finial states it re- 
duces to t + = and t~ = — s. 

If the quantum-statistical effects arc neglected then 
f L a q {Ea)f e b q {E b ) = f e N q {E q ). The integration over E a can 
be easily performed in this case and, combined with the 
A~ 2 (1 ; x a , x b ) prefactor, gives In the same approx- 
imation the last line of (|A6|) does not depend on the 
distribution functions and gives so-called 'reduced cross 



section' 



2tt 



1 f dipaj f dt„ 



Using Eq. (|A7[) and integrating over E q we recover the 
usual expression for the scattering reaction density: 



OO _ 



To take the quantum-statistical effects into account we 
need to express energies of initial and final states in 
terms of the integration variables. By energy conserva- 
tion E b = E q — E a . Therefore E a and E b as well as the 
related momenta \p a \ and \p b \ are completely fixed by the 
second and third integration variables. Next we consider 
the final states. By energy conservation Ej = E q — E t . 
It remains to express Ei in terms of the integration vari- 
ables. Let us choose the coordinates such that p a points 
along the a;-axis and p b lies in the xy-plane. Then the 
momentum transfer q also lies in the same plane. Its com- 
ponents are given by q = |g| (cos 9 aq , s'm9 aq , 0) where: 



cos 6 ag = 



2E a E q - s 



m\pa\ 



The components of pi can be written in the form pi = 
\pi\ (cos 9ai, sin 9 ai cos (p ai , sin ai sin tp ai ). Then the angle 
between the vectors pi and q is given by 

cos 9 qi = cos 9 aq cos 9 al + sin 9 aq sin 9 ai cos tp ai . 

Using energy-momentum conservation we can express Ei 
in terms of this angle and the integration variables: 



1 s 

Ei = 9 -1 ■ 2 fl [^i 1 + X i ~ X j) 

2 s + q z sin 9 q i 

+ \q\ cos9 q i(\{l, Xi, Xj) — 4xi<f / s sin 2 ( 



(A8) 



Note that for 4m 2 (f > X(s, m 2 , m 2 ) the difference under 
the square root in Eq. (|A8|) can become negative for some 
angles. This means that such scattering angles are for- 
bidden kinematically and should not be integrated over. 

Since Eq. (|A8[) implicitly depends on 9 a i it is conve- 
nient to use this angle as an integration variable instead 
of t. The integration measure in Eq. (|A6[) is then modi- 
fied according to dt 2\p a \\Pi \ sin# a i d9 a i. 

To calculate the scattering amplitudes we need the 
three Mandelstam variables, s is an integration variable. 
t is given by 

t = to 2 + to 2 - 2E a E l + 2\p a \\pi\ cos9 ai . 

The remaining one, u, can be inferred from the Mandel- 
stam relation s + t + u = Y^, . 
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Appendix B: Generalized optical theorem and 
thermal cutting rules 

The generalized optical theorem is a consequence of the 
unitarity of the S-matrix and can be seen as a consistency 
condition for the amplitudes to ensure conservation of 
probability. It can also be seen as a consequence of the 
Cutkosky cutting rules (79l - [8ll ] for the computation of 
the discontinuities of Feynman diagrams. As such it can 
be applied to unstable particles at any given order of 
perturbation theory. We may write it as 

-i [Ma^b({h}, {Pi}) - M* b ^ a ({pi}, {h}) 

= E(n/ ^W)*f(5>-J>) 

i \ il / j k 

x M„^({Ma?i,})>tU({Pi}> {».})• ( B1 ) 

The amplitudes Ai a ^b include all contributing diagrams 
(at a given order of perturbation theory) and the sum 
on the right-hand side is over all possible real intermedi- 
ate states i which contribute to Ai a ^b- The generalized 
optical theorem can be exploited to see explicitly why 
the RIS subtraction works. To this end, we apply it to 



V = w + w 



FIG. 18: One-loop contributions to £<j> —> £<f> scattering at 
0(h 4 ). The cuts through the internal lepton and Higgs lines 
yield the s x s, s x t and t x t-contributions to £<j> —¥ £<j>. The 
cuts through single internal Majorana lines yield the interfer- 
ence terms which contribute to the CP-violating parameter. 



the forward scattering processes £(j> —J- £(f> and £(f> — > £<f>, 
see Fig.[T51 above the energy thresholds s > m? Ni (the 
contribution by N% real intermediate states etc. can be 
addressed analogously) and s > {mg + m^) 2 . 

We include all possible graphs up to 0(h 4 ). Further- 
more, we sum Eq. (|B1[) over all internal degrees of free- 
dom of initial and final states and absorb these in the 
'effective amplitudes' defined in Sec.|ll] We get for the 
process involving particles 



2^{Yl M ^<t>} = J dll^(2n) 4 6{k+p-q)~ e 

dof 's 

+ J dn e q dnt(2n) i S(k+p-q-r)E i ^ 

+ J dU q dnt(2n fS(k +p-q- r)S t ^ . (B2) 

The amplitudes squared on the right-hand side contain 
the relevant graphs including the vertex and self-energy 





FIG. 19: One-loop contributions to £<j> — » £<$> scattering at 
0(h 4 ). Due to medium effects the particles in the loop can 
be on-shell. 



contributions, sec Fig.[TJ whose interference terms lead 
to CP-violation in the particle decay. In the same way 
one finds for £<J> <-> £<f>: 



23 { £ M r ^u) = J dU^(2nfS(k+p- q)E l 

dof 's 

+ J dll q dn$ (2nrS(k + p-q- r)3^_^ 

+ [ dU q dUf {2n)H{k +p-q- r)5~_^ . (B3) 



As a consequence of CPT we have ^M. e<p^ec/> = 
Y^-Mgz^gz. Therefore the difference of the left-hand 
sides of Eqs. (|B2|) and (|B3|) as well as that of the third 
terms on the right-hand-sides vanish. Subtracting the 
right-hand sides we obtain 



J dll e q dllf (2tt )*5(k +p-q-r) [E t ^ - S^J = 
dHf {2ir)H{k +p-q) [S^ Ni - E^ N J , (B4) 



as a requirement for a consistent approximation of the 
amplitudes compatible with unitarity and CPT. It is ob- 
vious that this cannot be satisfied if the scattering am- 
plitudes are in tree-level approximation meanwhile the 
decay amplitudes violate CP. Wc show in Scc. lIIII how 
Eq. (|B4|) can be satisfied by replacing the two-body scat- 
tering amplitudes 3 — > 3'. The solution for 3' amounts 
to subtracting the real intermediate state contributions 
from 3. In order to obtain an equivalent result for the 
Higgs decay at high temperature (i.e. for > me + AIi) 
we need to consider different processes since the ampli- 
tude lA^j^jy. cannot be obtained as cut of the graphs 
in Fig.1181 One could try to draw and to cut graphs for 
£Ni — > £Ni scattering, but the obtained cuts are lepton 
number conserving and would drop out in the difference 
Yl = Yi — Yg. Instead the relevant contributions ap- 
pear as 'thermal cuts' of the f-channel contributions to 
£4> — > £(f> depicted in Fig.[l]3 Since they exist only at 
finite density we need to use finite temperature 'circling 
rules', which can be seen as a generalization of Eq. (jBljl . 
to compute the imaginary part of causal n-point func- 
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FIG. 20: The thermal width can be obtained using (causal) 
finite temperature cutting rules. 




(a) 





(d) 




FIG. 21: Circlings contributing to the self-energy (a),(b) and 
vertex (c-e) CP-violating parameter for Majorana neutrino 
decay and for Higgs decay (f),(g),(c-e). The contributions 
by graphs (b),(f) vanish since Ni and Nj cannot be on-shell 
simultaneously for i 7^ j. Contributions by graphs (c-e) are 
suppressed if the cut is through an internal Majorana neutrino 
line. 



smallest time is always circled. We can use this equation 
together with the circling rules given in [22I . l82j to com- 
pute the imaginary parts of the graphs in Fig. [19] and the 
CPT-conjugated process £ cj> — > £<fi. Taking the difference 
of both we obtain, similar to Eq. (|B4p : 

J dill dUf^rSik +p-q-r) [S^ - S^J = 
dIi^)H{k +p-q) [E m ^ - S^J . 



We can also use Eq. (|B5|) to compute the thermal widths 
which cutoff the s- and t-channel resonances by cutting 
the self-energy graphs as shown in Fig. 1201 Furthermore, 
the thermal CP- violating parameters can be obtained us- 
ing thermal cutting rules, see e.g. Fig.1211 and [22j. Al- 
together, the concept of RIS subtraction can be general- 
ized to include quantum-statistical effects using thermal 
quantum field theory in the real time formalism and a 
complete set of reaction densities can be computed. Note 
however that inconsistencies are inherent out of equilib- 
rium and arise e.g. at higher order in the expansion per- 
formed in Sec. lIVl Similar computations where performed 
in [3TI . in the imaginary time formalism of thermal 
quantum field theory. 



tions in the real time formalism [821 ] : 

23{i !F^i^(xi,... 1 x a ,...,x n ;Zj S )^ = (B5) 

not all 

T X^j* _1 ' F >( a; l'---'£a'---' a; »; Z j) _ 

Xi Zj 

—i P< (x\ , . . . , , . . . , x n , Zj ) j> , 

where 'not all' means that not all should be circled at 
the same time. It was shown in [22( that causal n-point 
functions are the ones relevant for the computation of 
CP-violating parameters. The vertex x a with largest or 



Appendix C: Rate equations 

In this appendix we present some detailed intermediate 
steps in the derivation of rate equations with quantum- 
statistical terms. To obtain these we need to assume that 
the system is close to thermal equilibrium. Under certain 
conditions it is then possible to reduce Boltzmann-likc 
equations to a set of rate equations for the systems evo- 
lution. Assuming that the detailed conditions given in 
Sec. lIVl arc fulfilled and using Eq. ([3"T]) . we get for the in- 
tegrated Boltzmann equations with quantum-statistical 
terms: 



■pk;q _~ 



pk-q 



(CI) 



qpk 

N, 



%{2i:f S(q-p- k){±E N ^(l - /,)(! + /,) [ ^ fN ' eq - ( 



^1 +1*4, 



5*^(1 - mi + i ^fj^q - (e-^ - 1)]}(1 - f Ni ) 



This expression is still exact with respect to deviations 
from equilibrium and, taking into account results before 
Eq. f|32|) . it has an obvious expansion in (J./T. We can see 



that there will be contributions proportional to A.f N , = 
(f N , — /jy 9 ), to [ijT and proportional to Af N ji/T at lin- 
ear order. The coefficients of these terms were introduced 
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in Sec. lIVI (Eqs. (|3"3"|) and (JM])) and dubbed reaction den- 
sities. We find for the 0(fi/T) expansion of each of the 
two terms in Eq. (|C1[) . i.e. for the general collision term 
ofadoiV in Eq. dSJ): 



dr\ qpk " 



Nab^ab^N^ab^N 



P. 



pk;q 



AjW 
~1n 



"AT /at 

6 

1 N JN 1 

N 



■q~\%a7fJa + SbTfJb JIN 



Ma + Mb / S ab<->iV ( X _ £ /iv) 



T 



"AT 



(1-TO) 7 ^ } ' 



where we defined the decay reaction densities as 

(X 7a ?) = | dlilt b (27t) 4 5(q - fc - p)XE N f%U , 

(X 7A T) =(JC(1 - /^M) , 

with / o6 = (1 - Ca/a 9 - For a general process 

aN -s-> 6 (where we allow again for deviations from equi- 
librium in f N ) we find in the same way: 



dU qpk " 

U,±± AT .-. W — ' 



Nab^aN++b- r aN<->b 



pk;q 



^aN^b J N _,D 



A/, 



—b 



N 



Ma c f eq\ , <■ Mfc req\ D 

— jrli-Sa/a )+& jT-'fc )lb 



b J N 

Ma — Mb / ^ a N^>b (1 — ^ f n) 



T 



—iff 



(C3) 



where 



(*7 



pqfc 
oAfb 



(27r) 4 5(p + g-fc)XS h /^/ a& . 



(i T r)E(x(i-/;') 7 °), 

with, now, / a b = (/q 9 + ^w/fc 9 )- Note that these re- 
action densities will tend to zero in the zero tempera- 
ture limit because of their dependence on the distribu- 
tion functions. Using the general result (|C2|) we obtain 
the contributions to Eq. ((Clj) . and therefore to Eq. 



proportional to A/ 



Ni ■ 



sH 



2 A/ Ni \ 

contributions proportional to A.f N . 
dY Li 



req "7AT, 



lit 

T ' 



M£ 

T 



dz 
±3 



JV^ T 



( c <f>efl q 



/f) 



A/, 



or just proportional to [iijT: 
dY L , 



sH 



dz 



JV 4 ?±^ + "Ni^±U (1 - /jV 4 ) 



q\'Ni 



"JV, V 1 /AM 

This leads immediately to the result in Eq. (|36[) . Simi- 
larly, we get using Eq. (|C3|) for the contributions (|124l) 
by (anti-)Higgs decay to the rate equations: 



dYL 



z dz 

dYL, 



sH 



iz 



Af fi£ 



'-'tji^tiNi "F ^4>->-eNi A/at; 

5 J-eq 70 

"0 -'JVi 



M« / ^<f>^£Ni "F ^4>->eNi 



s% 



T 

dYL, 



(C^// - (1 - //))-7eg-70 



In addition we get for the 'extra' terms in Eqs. (|23|) and 



D. extra 



oc (1 ± 1)( "^' (1 - 7^)7^ > + (1 T D^d + c„)( "^' (1 - /Wl )7^ 



+ (1T1) 



[M_l ^ATj^l^ ~ ^£<p^±Nj 



^Ni 



(C4) 



This means that for dYn/dz only the first term and for To reduce the two-body scattering terms in Eq. ((771) 
dYjqJdz the last two terms contribute. The terms pro- we can use the same methods to find for the general ex- 
portional to the difference of Sjy.^^ and lead to 

contributions proportional to the CP-violating parame- 
ter. 
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pression, 



equation: 



JdIL k a ^(2^S(k + p- q -r)E ab ^ ij 

x [(i - r/ja - fhWj ± (i - efdo- - efflah] > 

the Q{n/T) expansion: 



dll 



~'ab<-+ij 



:(l-r/ )(l-^/6)/^[(l±l)±(e 
>c + M d ±(M +J^ (7 ^ ) + (1±1)<7 ^ 



-I) 



T 



(l±l)<(^/^+ 6 ^/ ( 



)>!, 



T 



T 



f e 

SI rp J % 



T 



where the two-body scattering reaction density was in- 
troduced in Eq. (j3"5)) . The last expression results after 
first order expansion in ji/T. 



Appendix D: 2PI effective action and self-energies 

In this appendix we derive one-loop contribution to 
the self-energy of the Majorana field as well as one and 
two- loop contributions to the self-energy of leptons. 

1. 2PI effective action 

The 2PI effective action is defined as a functional of 
the one- and two-point functions consisting of an infinite 
sum of all 2PI vacuum diagrams [64j]. In practice, its 
expansion can be characterized in terms of the number 
of loops appearing in each diagram: 

zr 2PI [s, A] = J2 ir&p, y, a] . 



(2) (3) 

The two lowest order contributions, il^pj and il^pj, rel- 
evant for leptogenesis are shown in Fig. [5] Their contri- 
butions to the 2PI action read 



iT 



(2) 



2 PI 



J d 4 u d 4 w r Tr[hP R y(u,w)P L /i f 

x S(w,u)eA*(u,w)e] , (Dla) 

* r 2Pi = o / diu dAyj A d ^ Tr[hP R ,y(u,w) 
2 Jc 

x CP R h T S T {r),w)eA*(w,Oeh*P L C 
x y(r),Z)P L rfS(£,u) eA*(u,?7)e] . 



In Eq. (jDip the trace is taken over flavor, Dirac and 
SU(2)l indices whereas the transposition only acts in 
flavor and Dirac space. 

2. Lepton self-energies 

By functional differentiation of the 2PI effective action 
with respect to the two-point function we obtain the cor- 
responding self-energy which enters the Schwinger-Dyson 



J a/3 



{x,y) = -i 



sr^[s] 



Here, flavor indices are shown explicitly whereas the 
SU(2)l and Dirac structure is embodied implicitly in 
matrix notation. The resulting self-energy is given by a 
combination of the Majorana, lepton and Higgs propaga- 
tors. Since we do not consider the flavor effects and the 
early Universe was in an SU (2) ^-symmetric state: 



A ab( x ,v) = $abA(x,y) 



5 ab S a0 S(x,y) 



Since e 4 = — e 2 = 1 the lepton self-energy also becomes 
diagonal: T,"^ = S ab T, a ^ . Furthermore, in the unflavored 
approximation it is convenient to sum over lepton fla- 
vors: E = ^ Q S aa . Then the one- and two-loop order 
contributions to the lepton self-energy, see Fig. 03 take 
the form: 



£ (1) (a;,y) = -(^h^Pn^ (x,y)P L A(y,x) , (D3a) 
& 2) (x,y) = -{rfh) l3 (rfh) lk f d 4 wd 4 V Pn^ jk (x,w)C 



x P R S T (r),w)P L C^ H (r),y)P L A(y,w)A(r),x). (D3b) 

Eventually, it is the Wightman components that we are 
interested in since they enter the gain- and loss terms on 
the right-hand side of Eq. (|55p. Therefore, we insert the 
usual decomposition of the propagators G £ {A, 3,^} 
into the spectral and statistical parts, Eq. (|4T|) . into the 
self-energies (|D3a[l and (|D3b[) . A formal decomposition 
of the self-energy in analogy to P?)) allows us to identify 
its spectral and statistical part and define the Wightman 
components in coordinate space as = S F =p ^E p . For 
the one-loop self-energy (ID3ap they read 



T^\x,y) = -{rfh) Tl P R ,^{x,y)P L AAy,x) 



(D4) 



In the case of the two-loop contribution, Eq. (|D3b|) . the 
computation becomes slightly elaborate. The complica- 
tion is due to the appearance of 32 different terms af- 
ter inserting the decomposition (|47]) for each of the five 
propagators into Eq. (|D3bp as well as due to the two 
remaining integrations over the internal space-time ar- 
guments w and 7/. The decomposition makes the path- 
ordering explicit and allows us to convert the integra- 
tion along the CTP into an integration along the positive 
branch. The 32 terms contain different combinations of 
the sign-functions. These can be rewritten by using rela- 
tions given in Appendix C of [2Q] • After some simple but 
lengthy algebra we obtain for the Wightman components: 
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T,f{x,y) = {h)h) l3 {h)h) lk J dojj d v 

x [P R .y^{x,u)CP R S T F {i h Lo)P L C.y^{i 1 ,y)P L A^(y,u)A A (r,,x) 
+ P R y 1 F k (x 1 Lo)CP R S T R (i 1 ,u)P L C,y^{i 1 ,y)P L A $ (y,uj)A A (ri,x) 
+ P R y R k {x 1 u J )CP R S T A {i 1 ,u)P L C^{i 1 ,y)P L A^(y,u;)A F (r],x) 
+ P R ^l k (x, u)CP R g£{ri, u)P L Cf]£{Ti, y)P L A F (y, uj)A $ (r), x) 
+ P R y^{x 1 Lo)CP R S T A (i 1 ,u)P L C,^{i 1 ,y)P L A R (y,uj)A $ (ri,x) 
+ P R y^ k (x,u J )CP R S];(r l ,u J )P L C^(r 1 ,y)P L A R (y,w)A^(r),x) 
+ P R y R k {x 1 u J )CP R Sl{i 1 ,u)P L C^{i 1 ,y)P L A^{y,uj)A^{ V ,x) 

+ P R yf(x 7 u J )CP R Sl( V ,u J )P L Cy^(r 1 ,y)P L A R (y,oj)A A (r,,x)] . (D5) 



Specific approximations will allow us to interpret both 
expressions, Eq. (|D4[) and Eq. (|D5|) . as describing decay, 
inverse decay and scattering processes of quasiparticlcs 
in the medium. 

The expressions for the one- and two-loop self-energies 
given by Eqs. (|D4[) and (|D5|) depend explicitly on two 
coordinates in four dimensional space-time. However, 
the self-energies which govern the gain- and loss term 
on the right-hand side of Eq. (|55|) are expressed in terms 
of phase space coordinates. Let us therefore exchange 
the pair of space-time arguments (x, y) for an equivalent 
set of center and relative coordinates, X = (x + y)/2 and 
s = x — y. In contrast to thermal equilibrium, the out 
of equilibrium propagators depend not only on the rel- 
ative coordinate s but also on the center coordinate X. 
Performing a so-called Wigner transformation [56j , i.e. a 
Fourier transformation with respect to the relative coor- 
dinate s, we can trade the latter for a momentum space 
variable: 

G F (X uv ,p)= /"d 4 Sw e* s -G F (X w , STO ), (D6a) 

G p (X uv ,p) = -i J d i s uv e ipSuv G p {X uv ,s uv ) , (D6b) 

where we have used X uv = (u + v)/2 and s uv = u — v, 
(u,v) <E {x,y,w,i]} according to the various combina- 
tions appearing in Eqs. (|D4j) and (|D5|) . Note that the fac- 
tor — i in the definition (|D6b[) is conventional and makes 
the Wigner transform of the spectral propagator a hcrmi- 
tian matrix. Definitions of the Wigner transforms of the 
advanced and retarded propagators coincide with that 
for the statistical propagator. 

The Wigner transform of Eq. (|E>4[) is obtained straight- 
forwardly: 

z^ ) (t,p) = -(rfh) ji J dni dnfanfSiq-k-p) 
xP R y^(t,q)P L AAt,k). 



Note that motivated by the homogeneity and isotropy of 
the early Universe we only indicate time-dependence of 
the propagators and self-energy, t = X® y . To obtain the 
Wigner transform of Eq. (|D5[) we will use an additional 
approximation: each of the Wigner transforms of the 
propagators we replace by G(X uv ,p) —> G(X,p). This 
means that we neglect the variations of X uv from the 
center coordinate X = X xy at which the self-energy is 
evaluated. Technically, it corresponds to a gradient ex- 
pansion to lowest order and therefore disregards all mem- 
ory effects. This can be compared to the 'Stofizahlansatz' 
within the usual approach to the Boltzmann equation. It 
is convenient to represent the resulting expression as a 
sum of three terms: 

E^E^ + E^+E^. (D7) 

The first term on the right-hand side of Eq. (|D7[) corre- 
sponds to the Wigner transform of lines one to three and 
four to six in Eq. (|D5jl : 

^\t,p) = J dl^ dlii (2tt) 4 % -q + k) (D8) 

x [(h^h) ij (h''h) lk A jk (t, q, k)P L Cy^(t, q)P L A^(t, k) 

+ (h^h) ji (h^h) k iP R y^(t,q)CP R V kj (t,q,p)A^(t,k)] . 

where we have introduced two functions containing loop 
corrections: 

A jk (t,q,k) = J dUt ± dUt 2 dUi 3 

x (2tt) 4 5((7 + fei + fc 2 ) (27r) 4 (5(fc + k 2 - k 3 ) 
x [P R y 3 R {t, -k z )CP R Sl{t, k 2 )A A (t, h) 

+ p R y F k {t, -k z )cp R s T R {t, k 2 )A A (t, h) 

+ P R yi k (t,-k 3 )CP R Sl{t,k 2 )A F {t,k 1 )] , 

and V kj {t,q,k) = P A\-(t, q, k) P. As we will see, (|D8|) 
describes CP-violating decay of the heavy Majorana neu- 
trino. The second term on the right-hand side of Eq. (|D7|) 
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is given by the Wigner transform of the seventh line of 
Eq. (|D5|) and describes lepton number violating scatter- 
ing processes: 

Zf 2 \t,p) = J dll 4 2 dUi i dn k2 (2n) i 5(p + h -P2- k 2 ) 

x (h^h) ij (h^h) lk [P R yi k (t,p 2 + k 2 )CP R Sl{t, -p 2 )P L 
x Cy^{t,p 2 -k 1 )P L A^(t,k 1 )A^(t,-k 2 )] . 

Finally the last term in Eq. (|D7[) corresponds to the last 
line of Eq. (|D5]l . 

E^ 3) (*,p) = / dn p2 dTl kl dTl k2 (2ir) 4 6(p + q 1 -p 2 -q 2 ) 

x (h^h) ij (h^h) lk [P R y^(t,-q 1 )CP R Sl{t,p 2 )P L C 
X S>g(t, q 2 )P L A A {t, -q 2 - p 2 )A R (t, q x - p 2 )} , 

and can be identified with lepton number conserving pro- 
cesses which do not contribute to generation of the lepton 
asymmetry. 



3. Majorana self-energy 

Differentiating Eq. (|Dla[) with respect to the two-point 
function of the Majorana neutrino and using definitions 
of the CP-conjugate two-point functions we obtain for 
the Majorana self-energy 

Ilv{x,y)=-g w [(rfh) ij P L S{x,y)P R A(x,y) (D9) 
+ (rfh)* j P R PS(x,y)PP L A(x,y)] , 

where we have assumed the SU(2) l symmetry of the 
medium and neglected flavor effects. The factor g w — 2 in 
Eq. (|D9|) comes from the summation over the SU(2)l in- 
dices. The CP conjugate self-energy differs from Eq. (|D9|) 
only in the propagators replaced by their CP conjugate 
counterparts and the couplings replaced by their complex 
conjugates. 

From Eq. (|D9j) we can read off the Wightman compo- 
nents of the self-energy: 

n |(z,y) = -9w[(rfh)ijP L S^(x,y)P R A^(x,y) 
+ (rfh^PRPS^y) PP L A^(x,y)} . 

Its CP-conjugate can be obtained by complex conjugat- 
ing the couplings and replacing the two-point functions 
by their CP- conjugates. To calculate amplitudes of the 
scattering processes we will need its Wigner transform: 

n|(<, q) = - g w J dni dn*(2ir) 4 s( q - P -k) 

x [(rfh) lJ P L S^(t,p)P R A^t,p) 

+ (rfh) ji P R PS>(t,p)PP L A>(t,k)] . (D10) 



From Eq. (|D10|) we can deduce the Wigner transform of 
the corresponding spectral self-energy: 

^(t,q) =- ^-[(h^h) lJ P L U p (t,q)P R 
+ (rfhy ij P R TL p (t,q)P L ], 
where we have introduced 

IL p {t,q) = 167r/ dU 4 k dU 4 {2Tr) 4 d(q-k-p) 

x [A F (t,k)S p (t,p) + A p {t,k)S F (t,p)] , 
U p (t, q) = WttJ dU 4 dll 4 (2ir) 4 6(q -k-p) 

x P[A F (t, k)S p (t, P ) + A p (t, k)S F (t, P )] P , 

and q = — q). In the quasiparticle approximation the 
Wigner transforms of the two-point functions of leptons 
and the Higgs are given by Eqs. (|56p ([5T)) and Eqs. ([M]) 
(|65p respectively. In a CP-symmetric medium, which 
the early Universe was to a very good approximation, 
A p (t,k) = A p (t,k) and A F (t,k) = A F (t,k). The ho- 
mogeneity and isotropy of the early Universe further- 
more imply, that there is no dependence on the momen- 
tum direction and the spatial central coordinate so that 
A p (t,k) = A p (t,k) and A F (t,k) = A F (t,k). Just like 
for scalars, in a CP-symmetric medium the fermion two- 
point functions are related by S p (t,p) — S p (t,p) and 
Sp(t'P) = S F {t,p). As for the p — > p transformation, 
the terms in the lepton propagators which carry spinor 
structure are not invariant under it: 

P L fP R = P L PfPP R = PPiifPhP ■ (Dll) 

bmce P 2 = 1, Eq. ([DTT|) implies that in a homogeneous, 
isotropic and CP-symmetric medium n and IT are left 
and right projections of the same 'vector' integral L p : 

n P J = -^W h )n p L + (^Wh] v 

Explicit form of L p depends on the kinematic regime and 
is presented below. 

To evaluate the decay amplitudes as well as amplitudes 
of the s-channel scattering processes we need to evaluate 
it for positive q 2 = M 2 and q 2 — s. If q° is also positive 
then L p takes the form: 

L p (t,q) = 16tt J dUtdU e p {2n) 4 S(q-k-p)f 

x [l + U(t,k)-ftQt,p)] . (D12) 

To obtain Eq. (|D12[) we have used the quasiparticle ap- 
proximation for the two-point functions and integrated 
over the zeroth components of the momenta. Integrating 
out the energy-momentum conserving delta-function we 
obtain for its Lorentz components: 

2P 

L° p = —h (yo,y) , 
m y 2 

- \{vl - y 2 )(! + x ? - x 4>) i q (vo,y)] i 



3G 



where j/o = <7° /T and y = \q\/T. The integral functions 
/„ arc defined by 

z + 

where, in complete analogy with Eq. (|A2[) . the integration 
limits are given by 

z ± = \[y Q (l + x t - xt) ±y\i(l,Xi,x (t> )] . (D14) 

For positive q 2 and negative qo the components of L p are 
related to the ones above by L° p {— q°, q) = L^(q°, q) and 

L p {-q°,q) = -L p (q°,q) respectively. 

To evaluate amplitudes of the t- and u-channel pro- 
cesses we also need to calculate L p for negative square of 
the momentum transfer. In this case momentum-energy 
conservation ensures that p° and k° cannot be positive 
or negative simultaneously. If they have different signs 
then, assuming homogeneity and isotropy of the medium 
and using relations f| 102[> . we find 

L p (t,q) = 1671-/ dXlf <m l p (2Ttff 

x{6(q + k-p)[f^t,k) + fi(t,p)] 
+S(q -k+ p)[U(t, k) + f- t (t,p)] } . (D15) 

Eq. (|D15|) implies that for negative square of the momen- 
tum transfer L p vanishes in vacuum. Although it is in 
principle possible to retain the thermal masses of leptons 
and the Higgs in the calculation, the resulting expres- 
sions are quite lengthy in this case. Neglecting the ther- 
mal masses we obtain for the Lorentz components of L p 
in this regime: 




where the integral functions are given by 

oo 

lt(y ,y)= f dzz n (—^— + ). 

n\yo,y/ j \e z + l e z Tva-iJ 

h(y±yo) 

Note that in this regime y > yo and therefore the lower 
integration limit is positive. 

To compute the scattering amplitude we need to cal- 
culate the product qL p . For the s-channel we find: 

qL p = q 2 (l + x t -x <j> )y~ 1 Io{yo,y), 

whereas the corresponding expression for the t- and it- 
channels reads 

qL p = q 2 y- 1 Y,Io(Vo,y)- (D17) 

± 



At low temperatures Eq. (|D17|) is exponentially small and 
vanishes in the vacuum limit. 

To analyze the Higgs decay we need to evaluate the 
spectral loop integral in a region of the phase space where 
the effective Higgs mass exceeds the sum of the Majorana 
and lepton masses. Using properties of the distribution 
functions under the p° — > —p° transformation we find 
after some algebra from Eq. (11181) : 

L p (t,q) = 16tt J <m k ^)H{q+p-k)f[f;" + ft cl ]. 

Just like in Eq. (|D12[) , the integration is over the (on- 
shell) momenta of the Higgs and lepton and the Majorana 
momentum serves as a constraint. Note that in this case 
we are interested only in the on-shcll Majorana momenta 
and therefore q 2 = M 2 . After integrating out the delta- 
function we obtain a result similar to Eq. (|D13[) : 

2T 

L P = — Ji (yo,y) , 
' y 

= T7^-[yo J i (yo,y) - 5(^0 -xt- i)J (yo,y)] , 
\q\ y 

where the integral function is defined as: 

The integration limits are given by an expression similar 
to Eq. (|D14jl but with 1 + xg — x^ replaced by x^, — X( — 1. 
When the effective Higgs mass approaches the kinematic 
limit, 77i0 = rri£ + M the upper integration limit ap- 
proaches the lower one, and the integral vanishes. 



Appendix E: Numerical parameters 

To perform the quantitative analysis we need to specify 
the Yukawa couplings. For simplicity we focus on the case 
of a very heavy third Majorana neutrino M3 3> M2 > 
Mi (MSM). In this limit the Yukawa couplings can be 
expressed in terms of the observed active neutrino masses 
and mixing angles and only one complex additional free 
parameter uj. In the Casas-Ibarra parameterization (83l — 
|86j the Yukawa couplings are given by 

Mi 

(fvh)u = -^-{m 2 |l - w 2 | + 777,3 |w 2 |) , 
v 

(h^h) 2 2 = ^r( m 3 U — uj2 \ + m 2 |w 2 |) , 

(tfh)i2 = Ml ^ 2 ( m2 wy/l - u 2 * - 777,3 Vl - w 2 ) , 



where v w 174 GeV is the Higgs vev and we have as- 
sumed normal hierarchy. In this case the physical neu- 
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FIG. 22: Washout parameter for iVi decay in the complex ui- 
plane. The values for the benchmark points are K\(w\) — 
47.20 and Ki(u)2) = 54.97. K\ has a minimum value of 
ifi(0) = 8.28. 



trino masses are given by 



m 2 
m 3 



(Am* o/ ) 3 k, 8.71 • lCT^ GcV 



(Am 2 



sol 



5.0-10" 



-li 



GeV. 



For illustration we choose the benchmark points uj = 
cxp(-0.01 • I) and w = cxp(-0.5 • I) denoted by BM1 
and BM2 in the plots. As masses of the right-handed 
neutrinos we choose Mi = 10 9 GeV and M2 = \/T0Mi 
respectively. This choice of parameters is such that ef- 
fects related to the resonant enhancement will be unim- 
portant but contributions from both heavy Majorana 
neutrinos can be relevant. Note however that there 



are lower bounds on the washout parameters Ki = 
Ti/ H\ T _ M = rhi/m^ (with 'equilibrium neutrino mass' 

m* = 167r 5 / 2 y^7« 2 /(3\/5A/p;)), see Figs.[22] and El 
The freeze-out of the asymmetry will therefore typically 
occur late (i.e. T <C Mi) which renders medium effects 
in general small. However for the qualitative issues dis- 
cussed in this paper our preference is to specify a consis- 
tent set of parameters for which we can discuss the gen- 
eration of the lepton asymmetry in terms of two heavy 
Majorana neutrinos. This is of course not a general re- 
striction for the employed techniques. 



K 2 



700 




FIG. 23: Washout parameter for N2 decay in the complex em- 
plane. The values for the benchmark points are Ka(u)i) = 
9.22 and ^2(^2) = 53.38. K2 has a minimum value of 
K 2 (-l) = K 2 (l) = 8.28173. 
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